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We use a combination of numerical density matrix renormalization group (DMRG) calculations 
and several analytical approaches to comprehensively study a simplified model for a spatially 
anisotropic spin- 1/2 triangular lattice Heisenberg antiferromagnet: the three-leg triangular spin 
tube (TST). The model is described by three Heisenberg chains, with exchange constant J, coupled 
antiferromagnetically with exchange constant J' along the diagonals of the ladder system, with peri- 
odic boundary conditions in the shorter direction. Here we determine the full phase diagram of this 
model as a function of both spatial anisotropy (between the isotropic and decoupled chain limits) 
and magnetic field. We find a rich phase diagram, which is remarkably dominated by quantum states 
- the phase corresponding to the classical ground state appears only in an exceedingly small region. 
Among the dominant phases generated by quantum effects are commensurate and incommensurate 
coplanar quasi-ordered states, which appear in the vicinity of the isotropic region for most fields, and 
in the high field region for most anisotropics. The coplanar states, while not classical ground states, 
can at least be understood semiclassically. Even more strikingly, the largest region of phase space 
is occupied by a spin density wave phase, which has incommensurate collinear correlations along 
the field. This phase has no semiclassical analog, and may be ascribed to enhanced one-dimensional 
fluctuations due to frustration. Cutting across the phase diagram is a magnetization plateau, with a 
gap to all excitations and "up up down" spin order, with a quantized magnetization equal to 1/3 of 
the saturation value. In the TST, this plateau extends almost but not quite to the decoupled chains 
limit. Most of the above features are expected to carry over to the two dimensional system, which 
we also discuss. At low field, a dimerized phase appears, which is particular to the one dimensional 
nature of the TST, and which can be understood from quantum Berry phase arguments. 



I. INTRODUCTION 

The nearest-neighbor spin- 1/2 Heisenberg antiferro- 
magnet on the triangular lattice is an archetypal model 
of frustrated quantum magnetism. While the isotropic 
model in zero field is rather well-understood and is known 
to order into a coplanar "120°" state, away from this 
limit the situation is less clear. Two deformations of the 
Hamiltonian are of particular physical and experimen- 
tal importance: the application of an external magnetic 
field and the introduction of spatial anisotropy into the 
exchange interactions. 

The spatial anisotropy is introduced by decomposing 
the lattice into chains with bonds of strength J, ar- 
ranged into a parallel array, with inter-chain interactions 
of strength J' (see Fig.[T]). Here, we define R=l — J' / J 
as the degree of anisotropy, and h measures the applied 
magnetic field. There have been many extensive stud- 
ies that consider these effects separately. However, a 
two-dimensional (2d) phase diagram, taking both effects 
together, remains to be understood. This problem is 
of considerable experimental interest. The application 
of a magnetic field is one of the few general means to 
tune quantum magnets in situ, and provides very impor- 
tant information on the quantum dynamics, as well as 
clues to the underlying spin Hamiltonian, which is often 
not well-known. Two materials whose behavior in mag- 
netic fields has been extensively studied are CS2CUCI4 
and Cs2CuBr4, which are known to be approximately de- 



scribed by the spatially anisotropic version of the model, 
with larger anisotropy in the chloride (R ~ 0.7) than the 
bromide (R ~ 0.3 — 0.5). Both materials exhibit a rich 
structure of multiple phases in applied magnetic fields, 
for which a theoretical view of the phase diagram would 
be quite helpful. 

The solution of the ground state of a fully two dimen- 
sional frustrated quantum spin model in a two-parameter 
phase space is quite ambitious. Here, we consider a 
somewhat simpler task, by concentrating on the prob- 
lem defined by the model confined to a cylinder with a 
circumference of three lattice spacings (i.e. making y 
periodic with period 3), which we refer to as a Triangu- 
lar Spin Tube, or TST (see Fig. |2j. By a combination 
of analytical approaches and extensive numerical simula- 
tions using the Density Matrix Renormalization GroupP 
(DMRG), we reveal a rich and complex phase diagram 
for the T ST, s hown in Fig. [3J We argue in the Discus- 
sion (Sec VIII) that much of this diagram translates to 
the fully 2d model. Whenever possible, we use a nomen- 
clature for the ground state phases which translates di- 
rectly to two dimensions, though there are, of course, 
differences due to the absence of spontaneously broken 
continuous symmetry in one dimension. 

Different parts of this phase diagram will be discussed 
in detail in the bulk of the paper, but we will highlight 
a few aspects here, where strong quantum features oc- 
cur. First, the isotropic line, R = 0, as a function of 
magnetic field has been considered many times in the 
two-dimensional limit. There, semi-classical methodd^ 
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FIG. 1: TST in (a) Cartesian and (b) sheared coordinates 
with intrachain interactions J and interchain interactions J' . 



predict the stabilization of both coplanar spin configura- 
tions by quantum fluctuations, and, most interestingly, 
a magnetization plateau, at which the magnetization of 
the system is fixed (at T = 0) at 1/3 of the saturation 
magnetization over a range of magnetic fields. We will 
refer to this state as the "1/3 plateau" throughout this 
text. On the plateau, the spins order into a collinear con- 
figuration. Stabilization of such a plateau is very much a 
quantum effect and is one of the more striking quantum 
features of the TST. The presence of the plateau has been 
confirmed for the two-dimensional spin- 1/2 Heisenberg 
model by exact diagonalization methods. Our DMRG 
study of the TST is also consistent with the semi-classical 
picture along the R = line. We directly confirm the two 
"coplanar" phases, and accurately locate the boundaries 
of the 1/3 plateau. 

Another regime of strong quantum fluctuations occurs 
when R is close to 1, where the system is composed of 
weakly coupled (strictly) one dimensional (Id) chains. 
There, an approach based on scaling and bosonization 
methods is possible, following Refs. I4|51 Those tech- 
niques (explained in this context in Sec.[V]) predict a spin 
density wave (SDW) state over a wide range of applied 
fields. In this SDW state, the dominant spin correlations 
are those of the Ising component parallel to the field, 
in sharp contrast to the classical behavior. Our DMRG 
simulations show that the SDW state dominates a re- 
markably broad region of the phase diagram, extending 
far beyond the decoupled line, R = 1. 

In two dimensions, the quasi-lD approach of Refs. I4|5l 
shows the existence of a (very narrow) 1/3 plateau arising 
out of the SDW phase, leading to the speculation that the 



plateau persists for all R in two dimensions. In the TST, 
we find that the plateau is also very robust, and persists 
almost, but not quite, to the ID limit. The suppression 
relative to two dimensions can be understood as a result 
of enhanced fluctuations due to the one-dimensionality 
of the TST. To check this, we have also carried out some 
DMRG studies of wider cylinders consisting of 6 and 9 
sites in the periodic direction. Our results appear con- 
sistent with the existence of a plateau for all R in two 
dimensions. 

The last quantum regime we discuss here is clearly spe- 
cific to the periodic boundary conditions imposed around 
the TST. This occurs at zero field, where for all values of 
i?, we observe a spontaneously dimerized ground state. 
The dimerization is most clearly observed in the entan- 
glement entropy, which shows a pronounced oscillatory 
behavior along the chain. We argue that this can be un- 
derstood as an effect of one-dimensional quantum fluc- 
tuations upon an underlying short-range spiral magneti- 
cally ordered state, somewhat similar to the formation of 
a Haldane gap in integer spin chains with collinear clas- 
sical states. The elementary excitations of the dimerized 
state are solitons, and we show how the behavior at small 
magnetization can be understood in terms of a dilute sys- 
tem of such solitons. 

The remainder of the paper is organized as follows. In 
Sec. [n] we introduce the model and then describe key 
technical aspects of our DMRG simulations, including 
the procedure to determine the phase boundaries using 
the second derivative of the ground state energy and en- 
tanglement entropy, and careful finite size scaling. In 
Sec. |III[ we review and compare the semi-classical pre- 
dictions to the DMRG results in the isotropic limit. Next, 



we discuss the high field region in Sec. [IV| In the vicin- 
ity of the saturation field, the problem can be modeled 
as a dilute system of spin-flip bosons. We compare an 
analysis of this limit, built upon an analytic solution 
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FIG. 3: (Color online) Phase diagram for the spatially 
anisotropic spin- 1/2 TST in a magnetic field. Here, we use 
the following abbreviations to label the various phases of the 
diagram: C — commensurate; IC — incommensurate; SDW — 
spin density wave. R = 1 — J' / J is the degree of anisotropy. 
The dashed lines indicate constant magnetization lines, where 
the upper, middle and lower ones are at M/M s = 5/6, 1/2 and 
1/6, respectively. 



of the Bethe-Salpeter equation, to the DMRG, and find 
a transition between coplanar and cone phases, and a 
commensurate-incommensurate transition. In Sec.[Vj we 
study the regime of weakly coupled chains, and in partic- 
ular discuss the spin density wave (SDW) state and show 
that the 1/3 plateau terminates in a Kosterlitz-Thouless 
transition around R ~ 0.7 ± 0.1 for the TST. We con- 
sider the low field region in Sec. |VII[ showing the per- 
sistent dimerization, the evidence for solitons at small 
magnetization, and the commensurate to incommensu- 
rate transition near R = 0. DMRG numerical results 
will be presented throughout these sections, presenting 
the important features used to identify each phase. Phys- 
ical quantities, like entanglement entropy, vector chiral- 
ity and the spin density profile will be shown for some 
representative large system size. Finally, we conclude in 



Sec. VIII with a summary and discuss some generaliza- 
tions of our results to larger spin and two-dimensional 
systems. 



II. MODEL AND DMRG METHOD 

A. Hamiltonian and notation 

The explicit Hamiltonian studied in this paper is writ- 
ten as 

H = y ] [J S x ,y ' &x+l,y + J' S Xj y • (S Xj y+i + Sa;_i j2/ -|_i)] 
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FIG. 4: (Color online) Second derivative of the ground state 
energy with respect to R, for different values of magnetization. 
These plots are used to locate the phase boundaries in Fig. [3] 



where x is the direction along the chains, and y is perpen- 
dicular to it, and h is the magnetic field. Importantly, 
we choose coordinates, as shown in Fig. [TJd, where the 
triangular lattice is "sheared" to embed it in a square 
one. This is convenient for the application of periodic 
boundary conditions in the TST. 

Many previous works on the anisotropic triangular lat- 
tice in two dimensions, including those by some of the 
authors^, use instead "cartesian" coordinates, as shown 
in in Fig. |TJl. Both for convenience in certain calcula- 
tions (especially in the quasi-one-dimensional limit), and 
to clarify the connection to this prior work, we give the 
relation between the sheared and cartesian coordinates 
here. In cartesian coordinates, we take the distance be- 
tween sites along the chains and the (normal) distance 
between chains to unity. Defining the cartesian coordi- 
nates as x, y, and r = (x, y), then 



x = x + y/2, 



y- 



(2) 



From this, we may also obtain the relationship between 
wave vectors in the two coordinate frames. We require 
q • r = q • r, which implies 



q x 



(3) 
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B. DMRG 

Throughout this paper, we rely extensively on DMRG 
simulations. For the present study, we kept up to m = 
3072 states in the DMRG block, performing more than 
24 sweeps to obtain fully converged results. In doing so, 
we find that our truncation error is of the order 10 -7 . We 
also take advantage of the cylindrical boundary condition 
to study large systems and to reduce finite-size effects 
for a more reliable extrapolation to the thermodynamic 
limit. In particular, in the regions above the 1/3 plateau, 
we find that observables have much better convergence, 
with a truncation error of the order 10 -9 . Even in the 
regions below the 1/3 plateau not close to the dimerized 
phase, we find reasonable convergence, with a slightly 
larger truncation error on the order of 10 _T . However, 
when we approach the dimerized phase near zero magne- 
tization, finite-size effects dominate: system sizes up to 
TV = 180 x 3 do not provide a reliable extrapolation to 
the thermodynamic limit. 

The phase boundaries in Fig. [3] were determined from 
the simulations. We describe the methodology for doing 
so here, leaving the characterization of the phases which 
occur for subsequent Sections. For the case of contin- 
uous transitions, it is common to calculate the second 
derivative of the ground state energy, d d ^§ •> as it will di- 
verge when the system undergoes a transition. For finite 
systems, however, one will observe a finite peak that in- 
creases with system size. For example, in Fig.[4Fa), where 
we show the second derivative for M/M s = 5/6 for vari- 
ous system sizes, we observe a transition at R = 0.6. This 
transition corresponds to the upper dashed line in Fig. |3J 
where there is a transition between an incommensurate 
planar and a cone phase. We use similar procedures to 
determine phase boundaries at other magnetizations, e.g. 
M/M s = 1/2, 1/6 in Figs. ||b,c) correspond to the mid- 
dle and lower dashed lines in Fig. [3j respectively. 

In addition to these divergent peaks, there are some 
other features (which are not phase transitions) due to 
finite size effects. For example, in Fig. [4ja) for M/M s = 
5/6, a broad peak near R = 0.8 actually decreases (and 
eventually goes to zero) in the thermodynamic limit. 
Therefore, we can confidently say that the cone phase 
dominates in the region R > 0.6, and that there is no 
transition at R = 0.8. Similarly, for Figs. [4jb,c), the 
fluctuations in the plots near R « 0.7,0.45, respectively, 
are finite size effects and vanish in the thermodynamic 
limit. 

Finally, we use the structure factor 

5AiAi (9) = ^E e " Jq ' (r " r ' ) ^r^>. (4) 
r,r' 

to determine the boundaries between the commensurate 
and incommensurate phases. For example, for small 
R, the transverse and longitudinal components of the 
structure factor peak at commensurate momenta Q = 
(47r/3, 27r/3) and (27r/3, 47r/3), respectively. This defines 



the "C planar" regions in Fig. [3j 

III. SEMI-CLASSICAL BEHAVIOR IN THE 
ISOTROPIC CASE 

A. Two-dimensional model 

The isotropic model, J' = J, has been extensively 
studied in two dimensions, and it is believed that a semi- 
classical description, with weak quantum fluctuations in- 
cluded via spin wave theory, is qualitatively correct in 
this case. We find that the semi-classical analysis largely 
carries over to the TST, with small modifications to allow 
for one-dimensional fluctuations. Therefore we review 
the established semi-classical results first. 

In the classical limit, where spins are described as 0(3) 
vectors, the isotropic problem is known to display an 
"accidental" degeneracy in a non-zero applied magnetic 
field 7 *. This can be seen from the fact that this model 
can be rewritten as 




where Sa = Si + S2 + S3 is the sum of the spins on a 
triangle, and the sum is over all triangles on the lattice. 
The ground state configuration is given by the constraint 

Sa-^z = 0. (6) 

At zero magnetization, this constraint is solved by plac- 
ing all spins in a plane, with the three spins in each tri- 
angle at 120° angles to one another in a three sublattice 
structure. A specific ground state is specified by three an- 
gles, e.g. two determining the plane of the spins and one 
determining the angle within the plane. All such states 
are related by 0(3) spin symmetry; so this is a symmetry- 
demanded degeneracy, a previous DMRG stud}S on the 
2d model also confirms the three sublattice structure. 

In a non-zero field, the ground states retain a three- 
sublattice structure, with three arbitrary angles remain- 
ing to determine the specific ground state. However, the 
presence of the field reduces the 0(3) symmetry to 0(2) 
(or U(l)), and only one of these angular degrees of free- 
dom is symmetry demanded. The remaining two angular 
degrees of freedom constitute an accidental degeneracy. 
Two simple states within the degenerate manifold are the 
coplanar and umbrella ones, shown in Fig. [5] 

As first shown by Chubukov and Golosov 2 , this ac- 
cidental degeneracy is lifted by quantum fluctuations. 
They showed by a 1/S spin wave expansion that the de- 
generacy is lifted in favor of the coplanar states. Ad- 
ditionally, they demonstrated the existence of the 1/3 
plateau, in which the spins adopt a 3 sublattice "up up 
down" structure. Away from the plateau, the copla- 
nar state retains a 3 sublattice structure with ordering 



5 



(a) 



(b) 



(c) 



cos Ox + sin Oy cos Ox + sin 6 



FIG. 5: Degenerate classical spin configurations in the 
isotropic limit. With the magnetic field taken in the z- 
direction, (a) shows the "V" configuration above the 1/3 
plateau, (b) depicts the "Y" phase below the 1/3 plateau 
while (c) shows the cone (or umbrella) state. 



wavevector Q = (4tt/3,0), or Q = (4tt/3, 2tt/3) 2 3 9 . Be- 
low the plateau, the 3 spins form a "Y" with one spin 
antiparallel to the field and two spins with equal positive 
projection to the field but at opposite angles from each 
other. This can be viewed as a deformation of the 120° 
state with spins in a plane containing the magnetic field. 
Here the spin configurations can be parametrized by 



(5+) = ae^sin(Q-r) 
(S z r ) = 6-ccos 2 (Q-r), 



(?) 



where is an arbitrary angle specifying the plane of the 
spins, while a, 6, c are constants dependent upon the field 
magnitude. Since Q r = 27r(2x + y)/3, we see from 
Eq. ^ that when 2x + y is a multiple of 3 the spin 
is aligned with the magnetic field, forming the vertical 
"leg" of the Y. Above the plateau one finds instead a 
"V" configuration, with two spins identical and the third 
chosen to give zero moment normal to z. In this case, we 
have 



ae i0 cos (Q • r) 
b — ccos 2 (Q • r 



(8) 



Note that the cosine in the first line of Eq. ^ never 
vanishes on lattice sites, so that spins are never parallel 
to the field in the V state. 



B. One dimension 

We will see that the semi-classical results summarized 
in the previous subsection for the two-dimensional case 
remain qualitatively correct, at least at short distances, 
in the TST. However, we must still account for the ef- 
fects of quantum fluctuations on long length scales, since 
the one dimensional system cannot break the U(l) spin- 
rotational symmetry about the field axis. Since the U(l) 
symmetry is unbroken in the plateau state, there are 
no essential effects of one-dimensional fluctuations there. 
However, they have qualitative effects in the Y a nd V 
phases, since (5+) = there, in contrast to Eqs. (7|8). 
Note that the modulation of (S*) is perfectly consistent 
with one-dimensionality, and is expected to persist di- 
rectly without qualitative modifications. 



To incorporate one-dimension al flu ctuations, we regard 
the semiclassical results in Eqs. ( 7|8 ) as defining the local 
spin ordering, with a fluctuating quantum phase 9(x,t) 
(r is imaginary time), that is, we make the replacement 

5 r + (r)^ae^' r) sin (Q • r) , 



(9) 



in the Y phase, and 

S+(f)^oe' 9(l; ' r » cos(Qt), 



(10) 

in the V phase. Note that these formulae are not in- 
variant under translations, reflecting the three-sublattice 
structure of the coplanar phases. This can also be seen 
from the oscillations in the (S*) expectation values. Even 
when one-dimensional fluctuations are taken into ac- 
count, translational symmetry is broken. This is still 
consistent with the Mermin- Wagner theorem, since the 
broken translational symmetry is discrete. Translating 
by one or two lattice spacings, one obtains two other 
symmetry related but distinct ground states. 

In both the Y and V phases, the field Q(x, r), represent- 
ing the "would-be" Goldstone mode of the spontaneously 
broken U(l) symmetry, is governed by the usual massless 
free relativistic boson action, 



S 6 



= J dxdT ^ 



{d x ef 



K 

2v 



(11) 



C. Comparison to DMRG 



We now turn to a comparison of the semi-classical pre- 
dictions, corrected as in the previous subsection for one- 
dimensional fluctuations, to the DMRG. 



1. Entanglement entropy 

The simplest comparison arises immediately from 
Eq. (11): the low energy physics is that of a single mass- 
less scalar field, which is a conformal field theory with 
central charge c = 1. This central charge can be directly 
measured using the entanglement entropy. 

According to conformal field theory^, in a one dimen- 
sional critical system with open boundary conditions and 
total length L, the von Neumann entanglement entropy 
associated to a region with length x and its complement 
of length L — x is given by 

L ■ m\ 



S(x,L) = -\n 
6 



sin 



(12) 



By plotting the entropy S(x,L) versus the reduced coor- 
dinate x' = ln[^ sin (^f )], we can directly extract c from 
the numerics. As shown in Fig. [6j we can indeed ob- 
tain c = 1 with high accuracy for both Y and V phases. 
For example, the obtained central charge c = 0.98 at 
M/M s = 1/6 in the Y phase below the plateau, and 
c = 0.97 at M/M s = 1/2 in the V phase above the 
plateau. Both are consistent with the theoretical pre- 
diction. 
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FIG. 6: Entanglement entropy in the isotropic limit, R — 
for system size iV^ = 120. Note that the reduced coordinate 
x = In [^sin(^)] is plotted on the x-axis. We show the 
von Neumann entanglement entropy for (a) M/M s — 1/6, 
the commensurate Y phase and (b) M/M s — 1/2, the com- 
mensurate V state. The solid line is a linear fit, where by 



Eq. (12), we can extract the central charge, c. 



FIG. 7: (Color online) S z profile in the isotropic limit, R — 
at (a) M/M s — 1/6, the commensurate Y phase, and (b) 
M/M s = 1/2, the commensurate V state. The square (red), 
diamond (green) and circle (blue) data points show the 3- 
sublattice structure of the isotropic case. The magnitude of 
S z does not decay because the discrete translational symme- 
try is spontaneously broken in these states. 



2. S z profile 

The modulation of (S*) predicted by the semi-classical 
theory in Eqs. ( |7|8[ ) can be directly compared to the 
DMRG results. This is shown in Figs. [T[HJ Note that 
a particular symmetry broken state is chosen in the sim- 
ulations, presumably due to pinning by the boundaries, 
which explicitly break translatio nal symmetry. The ori- 
gin of the coordinate r in Eqs. ( 7|8 ) must be appropri- 
ately chosen to match the chosen ground state. 



3. S correlations 

Due to quantum fluctuations of the phase 0, the single 
spin expectation value (S^~) = 0. Therefore, we must in- 
stead turn to correlation functions to detect the Y and V 
structure of the local ordering. Using Eq. ([9]), we obtain 



where 



(S+S~,) ~ a 2 sin (Q • r) sin (Q • r') (e 



,i{0(x)-0{x')) 



in the Y phase below the 1/3 plateau. A similar formula, 
with the sines replaced by cosines, describes the correla- 
tion function of the V phase above the plateau. The cor- 
relation function is evaluated with respect to Eq. (11), 
where a finite-size form, first derived in Ref. 11, is as 
follows 



[f(2x)f(2x')} 



\f(x-x>)f(x + x>)]*' 
2(L + 1) 



(15) 



■ sin 



Here ao is a cut-off dependent factor, which we can take 
to unity, absorbing the dependence in a in Eq. (13). The 



(e i(0(*)-0(s'))) = C V (X,X'), 



(14) 



function, /(#), originates from a quantum average over 
the normal modes of the bosonic field 6. One is now able 
to fit the DMRG measurement of the transverse spin- 
spin correlation function to Eqs. ( |13|14 ) to obtain the 
ordering wave vector and the additional fit parameter, 
rj. A comparison is plotted in Fig. [9j where we show the 
correlation function along each chain (i.e., y = 1,2,3) 
for R = and M/M s = 1/6, 1/2. The fitting in Fig. |9k 
yields a commensurate wave vector Q = (47r/3, 27T/3) 
and r] = 0.65 for M/M s = 1/6, which corresponds to 
the Y phase below the plateau. Above the plateau, in 
the V phase shown in Fig. [9}d, the ordering wave vector 
still shows commensurability, Q = (47r/3, 27r/3) with r] = 
0.43. One can show that in the thermodynamic limit, 
the correlation function in Eq. ([l4| reduces to a simple 
power-law relation cx \x — a/| _r7 , which is reflected by our 
data for distances \x — x'\ <C L/2. 



D. Behavior for small non-zero R 

If we perturb slightly away from the isotropic limit, i.e. 
< R <C 1, we expect the semi-classical picture to still 
hold. This has been analyzed in Refs. l3|9l Classically, the 
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FIG. 8: (Color online) S z profile for the commensurate V 
phase at M/M s = 1/2 and R = 0.1. We find that the wave 
vector remains commensurate, even for a non-zero, but small 
R. 




minimum energy spin configuration changes immediately 
when R > from a commensurate state to an incommen- 
surate one, with an ordering wavevector Q ^ (47r/3, 0) or 
Q =^ (47r/3, 27r/3). However, we expect that quantum 
fluctuations will stabilize the commensurate state for a 
range of anisotropics for a generic value of the magnetic 
field. The reason is that coplanar phases break discrete 
translational symmetries of the lattice. Since there are 
three equivalent ground states connected by translations, 
the symmetry breaking can be described by a Z3 order 
parameter. Specifically, the combination 



r _ qz 27ri(x+2y)/3 



(16) 



defines a Z3 order parameter with (£ r ) = \(\e and 
$ = 0, 27r/3, 47r/3 in the three distinct Z3 domains. To 
restore this discrete symmetry, a phase transition is re- 
quired. More specifically, there are topological exci- 
tations of the coplanar state which are domain walls, 
also called solitons, connecting different symmetry bro- 
ken states. There is a non-zero energy gap to create a 
domain wall in any phase with long-range Z3 order. For 
the Z3 order to be destroyed, solitons must proliferate in 
the ground state. Small changes of parameters, such as 
i?, cannot instantly lower the gap for the domain walls to 
zero, which implies stability of the phase for a range of 
R values. This is correct, at least, away from the excep- 
tional points where h = (where the symmetry breaking 
becomes continuous) and h = h sat (where the symmetry 
breaking vanishes). We will discuss the vicinity of these 
exceptional points in subsequent sections. 

In general, with increasing anisotropy, i?, we will en- 
counter a phase transition to an incommensurate phase, 
which corresponds to the proliferation of solitons and a 
vanishing of their gap. Beyond that point, (£ r ) becomes 
zero, and S z correlations peak at a wavevector other 
than Q = (47r/3, 2ty/3). This transition is discussed in 
Sec. OVeI 

A useful test for this phase is the measurement of the 
central charge via entanglement entropy. In the commen- 
surate regions, even for R > 0, we expect c = 1, while 



FIG. 9: (Color online) Transverse spin-spin correlations in the 
isotropic limit, R = in the (a) commensurate Y phase and 
(b) V phase for N x = 120. Data points are shown as (black) 
circles while the theoretical fit from Eq. ( 14 ) is shown as the 
(blue) line. 



incommensurate phases may have c > 1. We observe this 
effect in Fig. |6j which shows c = 1 in the commensurate 
state, whereas Fig. 11 shows c = 2 in the incommensu- 
rate state. In addition, we can check for commensurabil- 
ity using structure factor measurements, as discussed in 
Sec. HTBl 



E. Phenomenological analysis at low field 

We now address the region slightly away from R = 
and at low applied magnetic field. We begin the discus- 
sion from a 2d point of view, though it largely applies to 
the TST as well. Commensurate coplanar spin order is 
described by the order parameter d = ni — zri2, where ni, 
112 are mutually orthogonal vectors with identical norm 
spanning the plane of the spin order. Then, a spin at 
coordinate r can be written as 

S r = M + Re(de iQ " r ) = M + ni cos[Q • r] + n 2 sin[Q • r]. 

(17) 

Lattice translations transform d de * 27r / 3 ? while lat- 
tice inversion, r — » — r, results in complex conjugation, 
d — >• d*. The effective Ginzburg-Landau Hamiltonian de- 
scribing the coplanar state should remain invariant under 
these operations (see Ref. 12 for a closely related discus- 
sion). Then, 



-rd* • d + a |c^d| 2 + ai(d* • d) 2 + a 2 |d • d| 2 
+Xi^d*-d + X 2|h-d| 2 



1 



X 3[(h-d) 3 + (h-d*) 3 ]. 



(18) 



Here, at mean-field level, r > is required to obtain 
non-zero ni 5 2, and ao,i > 0, for stability in the ordered 
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phase. Furthermore, CL2 > energetically imposes the 
orthogonality condition ni • 112 = in zero field. To 
favor coplanar (rather than umbrella) spin structures in 
a finite magnetic field, requires X2 < 0. We may expect 
that X2 is a function of the anisotropy, being negative for 
the isotropic limit R = and changing sign to positive 
values for sufficiently large i?, where the order by disorder 
physics favoring coplanar states gives way to the classical 
energetic preference for umbrella states. Here we restrict 
ourselves to the small anisotropy regime, for which we 
expect X2 to remain negative. With the preference for 
coplanar states set by X2 < 0, for field oriented along i, 
the preferred configurations of d may be parametrized as 



d = \d\e iG [z + i(cos #x + sin <9y)] , 



(19) 



where describes the orientation of the plane of the spins, 
and the angle of the spins within that plane. With this 
form for d, we obtain the spin operators as 



s. 



M 



-\d\e~ 



I cos(Q • r - 
sin(Q • r • 



0), 
-6). 



(20) 



The last term in Eq. (18) describes the commensu- 



rate locking of the spin to the lattice by the finite mag- 
netic field. Using Eq. (19), it may be rewritten as a 
sine-Gordon term 



H ag =x 3 |C/i d cos[30]. 



(21) 



The sign, %3 > 0? is fixed by the condition that one of 
the three spins in a sublattice must be oriented opposite 
to the external field in the commensurate state. Thus, in 
the commensurate state, = tt in (19). 



Now we move away from the isotropic line to R > 0. 
Here 3-fold rotational symmetry is broken, which allows 
the introduction of an additional term, linear in deriva- 
tives, into the effective Hamiltonian: 



-^(d^d-d-cU*) 
= -b 1 \d\ 2 dJ. 



(22) 



Since this term must vanish at R = and be analytic, 
bi ~ R. This term competes with the sine-Gordon term 
in Eq. (21), with the commensurate state with constant 



favored at small R and destabilized at larger R. Thus 
the commensurate-incommensurate transition in two di- 
mensions can be described the Hamiltonian of the phase 

#c-ic = J d^idoidjf-hdj^xs^co^e}}. (23) 

Here, the coefficients with tildes, ao?^i?X3? are rescaled 
by unimportant factors, such as the amplitude \d\. 
The sine-Gordon model of the form in Eq. (23 ) appears 



in several guises in this paper, and is analyzed in Ap- 
pendix [Aj It encodes a commensurate-incommensurate 
transition (CIT) with increasing b\. This transition is 



mean-field like for d = 2, and we may apply the results 
of Appendix A 1 This gives a critical value for the CIT 
of 6i )Cr ~ ^doX3h 3 for the incommensurate state, which 
translates to 



ft, 



c-ic 



R 2 '\ 



(24) 



since b\ ~ R. This is roughly consistent with shape of 
the boundary in the lower left corner of Fig. [3] 

For the TST, the situation is complicated by one- 
dimensional fluctuations. At zero field, h = 0, we 
know that, in fact, the ground state is not a spiral but 
rather a dimerized phase. Hence, we cannot directly ap- 
ply the above analysis at the lowest fields. The dimer- 
ized phase is broken fairly rapidly by the field, and so, 
above some small critical field, we may expect to be able 
to use results of this type. Even so, we should really 
use results for the d = 1 case, where a non-mean-field 
anal ysis a pplies, as described in Appendix |A 2| Using 
Eq. ( A16| ), the critical value bi^ CY is suppressed by a 
factor of (x3^ 3 /«o) A3/(4_2A3) , so that the net result is 

3-Ao 

61, ' 



h*= 



and hence 



h, 



c-ic 



i? 3 



(25) 



Here, A3 is the scaling dimension of the cos 30 term. As- 
suming the commensurate phase is at all stable for small 
R implies A3 < 2, so that the cosine term is relevant 
in the isotropic case, R = 0. It is also bounded below 



by zero, so that the exponent in Eq. (25) varies between 



and 2/3. Once again, we caution that the expression 
must be taken with care, since it does not in fact apply 
at the lowest fields. 



IV. HIGH FIELD REGION 
A. Spin flip bosons 

In this section, we study the phase diagram near sat- 
uration, i.e. for applied fields sufficiently large that the 
magnetization is close to its maximum of 1/2 per site. At 
saturation, the ground state of the model is the trivial 
product state with all spins aligned in the direction se- 
lected by the field. For fields above the saturation field, 
this is the exact ground state, and the lowest excited 
states consist of single magnons, in which just one spin 
has been flipped relative to the saturated state. These 
magnons are bosons with S z = 1, and upon reducing 
the field to the saturation value, the minimum energy 
required to create a magnon vanishes. Below the satura- 
tion field, therefore, we can expect Bose-Einstein conden- 
sation (BEC) of these magnons. In the one-dimensional 
TST, strict BEC is not possible due to phase fluctuations, 
but these fluctuations are readily taken into account and 
a quasi-condensate description remains appropriate. 

To formalize the magnon BEC picture, one may trans- 
form the spin model to a bosonic one^"^, using the 
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equivalence of the spin 8 = 1/2 Hilbert space to that 
of hard-core bosons: 



5+ 



V r b r V r 
1 



(26) 
(27) 



where n r = 6j6 r is the boson occupation number, and 
one must project onto the space of no double boson occu- 
pancy, V Y |n r 0)(n r 0| + |n r = l)(n r = 1|. Eq. ([26 ) 
is equivalent to the Holstein-Primakoff bosonization for- 
mula, truncated to quadratic order in boson operators 
and taking s = 1/2, provided the no double occupancy 
constraint is imposed. The generalization to s > 1/2 will 
be briefly discussed later in Sec. |VIIIB l[ 

It is convenient to implement the no double occupancy 
constraint by first relaxing the constraint, adding an on- 
site interaction U to the Hamiltonian, and then realiz- 
ing the projection by taking the U oc limit. In this 
way we can proceed simply by rewriting the Heisenberg 
model using Eq. (26), forgetting the projection operators, 
i.e. taking V r — >> 1. We thereby obtain a boson Hamil- 
tonian with hopping terms (J), on-site energies (J, h), 
an on-site (U) and nearest-neighbor (J, J 7 ) interactions. 
Fourier transforming to diagonalize the quadratic terms, 
we find 



H= $>(k) 



2N 



where 



6(k) 

V(k) 



k,k',q 



Jrr 



J(k) - , 

ft-sat - h, where h sat = J (0) 
2(e(k) + U). 



(28) 



(29) 
(30) 
(31) 



Here, J(k) is the Fourier transform of the exchange in- 
teraction, \i is the bosonic chemical potential, and h sat is 
the saturation field. 



B. Effective field theory for dilute bosons 

For h > /i sat , the vacuum is an exact ground state of 
this Hamiltonian, i.e. 6k |0) = 0. Below the saturation 
field, a finite density of magnons is introduced into the 
system, and a BEC or quasi-BEC is expected. The phase 



of the system, and correspondingly the magnetic order 
(correlations) , is determined by the structure of this con- 
densate (or quasi-condensate). To determine this struc- 
ture, we construct an effective model. The lowest energy 
magnon excitations in the triangular lattice occur at non- 
zero momenta ±Q, which minimize the dispersior J 16 * 17 *. 
In our (sheared) coordinates, the dispersion relation is 

JxsT(k) = Jcosk x + J' [cos k y + cos(k y — k x )\. (32) 
In two dimensions, we can choose arbitrary k x and k y , 
and the minima occur at k = ±Q2d, with Q2d = 
(<32d,<22d/2), and 



Q2d = 2 arccos 



r 

2J 



(33) 



Note that in the conventional cartesian coordinates this 
wavevector is Q = (Q2d, 0). For the TST, we must quan- 
tize k y = 0, 27r/3, 47r/3. With this restricted choice of 
k y , the 2d wavevector Q2d cannot generally be achieved. 
Instead, we find that the minimum energy wavevector is 
k T ST = ±Qtst = ±(Qi d ,27r/3), with 



{Id 



7r + arctan 



V3J' 
2J-J' 



(34) 



The two wavevectors coincide when J = J' . 

In a low-energy description, the modes away from these 
two minima may be integrated out, leaving an effective 
theory in terms of two "flavors" of bosons, ipi and ^2, 
defined via 



6 k = ip : 



l,Q+k 



■ $ 



2,-Q+k 



•6k- 



(35) 



Here, 7/>i ?q (^2,q) is defined as a boson "centered" on the 
minimum energy momentum Q (— Q), with weight only 
for small \q\ < A, where A <C 2tt is a cut-off introduced 
by integrating out the modes away from the minima. 
The third operator 6k represents the high energy modes 
which remain uncondensed, and are integrated out. In 
two dimensions, Fourier transforming in q x ,q y back to 
real space leads to slowly varying continuum fields ^ a ( r )> 
where r is a two dimensional spatial coordinate. For the 
TST, we need to keep only the mode with minimum en- 
ergy q y: and so, we Fourier transform only in q x , which 
leads to a continuum field dependent only on the position 
along the chain, x. 

In this continuum limit, the boson fields are governed 
by an effective action of the form 



J d d TdT I 4(8 T - ^V 2 )^! + 4(d T - ^V 2 )^ 2 — (l(pi + P2 ) + { P \ + pj) + T 2Pl p 2 



(36) 



where p a = \ip a \ . We have written the action, Eq. (36), in a form which includes both the TST (d = 1) and two 
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dimensional (d = 2) cases. We expand to fourth order 
in \i/j a \ and to lowest order in derivatives, which is justi- 
fied near saturation due to the diluteness of the magnons. 
The quadratic terms in Eq. ( 36 ) can be readily extracted 



from the exact single-magnon dispersion, which is given 
in Eq. (29) (in general in two dimensions the quadratic 
term may have an anisotropic effective mass tensoi^, 
which is not explicitly shown in Eq. (36)). The quartic 



interaction terms are more subtle, because though the 
magnons may be assumed dilute, the lattice-scale inter- 
actions in Eq. (28) are not weak. Therefore the parame- 



ters Ti, T2 must be obtained from a more careful analysis, 
which we return to below. 



C. Order parameter structure 

Taking for the moment the T a as phenomenological 
parameters, we discuss the structure of the condensed or 
quasi-condensed phase. If p < 0, there are no bosons in 
the system, and the vacuum is the ground state. When 
H > 0, a finite density of bosons is present. Depending 
upon their interactions, different phases may resuliP^. To 
discuss the nature of these phases, a mean field analysis of 
Eq. (36) is sufficient. We comment on the modifications 



to the mean field results at the end of this subsection. 



In mean field theory, we simply minimize S in Eq. ( 36 ) 



for constant values of ip a . When p > and Ti < I^, 
then pi 7^ 0, p2 = or vice versa, which means that the 
magnons condense at one of the two minima: a single-Q 
condensate. Here, in minimizing the energy, one finds 
that p 1 = (pi) = p/T 1 and E/N = -/i 2 /(2ri). By 



taking ^1,2 
as follows 



o™l,2 

1 

St = 
SI = 



one can write the spin operator 



^ e *(Q-r+fll) 
\ ~ <Pl), 



(37) 
(38) 



where ip = \f{pi) in mean field theory. We see that 
the z-component of the spins is non-zero but constant in 
space, while the xy components rotate as one moves in 
space. Such a configuration is called a cone or umbrella 
phase, because the spins trace out a cone as one proceeds 
through the lattice, see Figure |5|c). 

When T2 < Ti, then pi = p2, which means that the 
bosons condense at both +Q and — Q. This is a double-Q 
condensate with density (p) = {pi) J r{p2) = /V(ri+r 2 ) in 
mean field theory. Here, the energy E/N = p 2 /(Ti+I^). 
Again, by letting -0 lj2 = J~p^e ldia and 61 ^ = ±6, 



5+ = 2^cos(Q 
1 



S* = 



-0 

4(p) cos 2 (Q • r + 9 



(39) 
(40) 



where x/j = y (p) in mean field theory. In this phase, the 
z-component of the spins is not constant, but the phase 
of is constant. This implies that the spins remain 



in a plane, i.e. this is a coplanar phase. Instead of a 
cone, the spins in this phase sweep out a "fan" - so this 
is sometimes called a fan state. 

How much of this survives beyond mean field theory? 
In general, the dependence of the density on chemical 
potential is affected by fluctuations. Note that in the 
original spin problem, this dependence gives the behav- 
ior of the magnetization versus field in the vicinity of 
saturation, as is seen from Eq. (26). As is well-knowrP^, 
the BEC transition at p = is a very simple example of a 
quantum critical point, whose upper critical dimension is 
d = 2. Thus in two dimensions, the deviations from mean 
field theory are minimal and consist just of logarithmic 
corrections. However, in d = 1 the corrections are much 
more significant, and the dependence of the density on 
chemical potential is quite different. 

In mean field theory, we see that there is a first or- 
der transition between the cone and fan states upon 
varying I\ — T 2 through zero. In fact, the location of 
this transition at Ti = T2 is correct and moreover, ex- 
act, beyond mean field theory. To see this, note that 
when Ti = T2 = T, the interaction terms may be 
rewritten as ^(pi + P2) 2 , which implies that the ac- 
tion has an enlarged SU(2) symmetry under rotations 
~- ^ ^ZpUapi/jp, where U is an arbitrary SU(2) ma- 
trix. This guarantees the degeneracy of the cone and fan 
states at this point, since one can be rotated into the 
other by such an SU (2) rotation, and therefore, fixes the 
location of the cone to coplanar transition. 

When Ti 7^ T 2 , the SU{2) symmetry of Eq. ([36]) 



is reduced to U(l)xU(l), corresponding to independent 
phase rotations of ^1 and ^2. As a consequence, there 
will be one gapless mode in the theory described by 
Eq. (36) for each bose field with non-zero amplitude, i.e. 



one in the cone state and two in the fan. The fluctuations 
of these gapless modes lead, in the one dimensional TST, 
to power-law correlations of the spin components trans- 
verse to the magnetic field, rather than the long range 
order (broken symmetry states) obtained in mean field. 

Physically, the overall U(l) symmetry under simulta- 
neous and equal rotations of both fields reflects conser- 
vation of S z , and is microscopically mandated by the 
Heisenberg model. The "orthogonal" symmetry under 
the rotation of the two boson fields by opposite phases 
is emergent, however. It is a consequence of the discrete 
translational symmetry of the lattice, and the (generi- 
cally) incommensurate nature of the wavevector Q. In 
general, this symmetry is broken by terms (which should 
be added to S in Eq. (36)) of the form 



(41) 



where naively q n = 2Q, but in fact we can take q n = 
2Q — K/n, where K is any reciprocal lattice (RL) vector, 
since r is a lattice coordinate. So henceforth we work with 



q n = mhiKeRLpQ - K/n], 



(42) 
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9 + 




k'-p 



k + q 



k' -q 



FIG. 10: Ladder approximation of Eq. (44). Here, k, k' are 



incoming momenta while k + q, k' — q are outgoing momenta. 



i.e. we choose K to minimize the magnitude of q n . When 
the wavevector Q is incommensurate and the magnitude 
of these terms are small, their oscillations average to zero 
over short distances, and they can thereby be neglected. 
However, if 2nQ is close to a reciprocal lattice vector, 
then q n is small and the corresponding w n term becomes 
slowly varying, and it can have effects that persist into 
the continuum theory. This occurs only if 2nQ is close to 
a reciprocal lattice vector and the amplitude of both ipi 
and ip2 is non-zero, i.e. within the coplanar or fan state. 
This leads to commensurate-incommensurate transitions, 
discussed in Sec. IIVEI 

In the cone state, such effects are not important. In 
this case we expect one gapless "Goldstone" mode (6\) 
and power-law transverse spin correlations. But actu- 
ally there is some hidden long range order. Note that 
in Eq. (37) we have (arbitrarily) chosen the minimum 
with pi ^ and p2 = 0, instead of the one with pi = 0, 
p2 7^ 0. In doing so, the system spontaneously breaks dis- 
crete symmetries. In particular, for the TST, this choice 
breaks both inversion symmetry and a "charge conjuga- 
tion" symmetry, the latter being the anti-unitary sym- 
metry of the Scrodinger equation under complex conju- 
gation of the wavefunction. Although the fluctuations of 
the phase 0i above will reduce the mean field magnetic 
order to quasi-long-range order in the TST, the discrete 
symmetry breaking is robust to one dimensional fluctu- 
ations. This symmetry breaking can be most directly 
sensed by the vector chirality^^U, 



(43) 



Replacing 5+ in Eq. (43) by the ansatz in Eq.(37), we 

2 

find V = ip sinQ, i.e. a non-zero and constant value in 
the cone state. The opposite sign would be obtained for 
the solution with pi = 0, p2 ^ 0, so this serves as an 
Ising-type order parameter for the cone state. 



D. Incommensurate planar to cone state transition 
at the saturation 



1. Bethe-Salpeter equation 



Now that we have described the phases of Eq. (36), 



we will briefly outline the methods to compute Ti,!^. 
When the external field is sufficiently close to the satu- 
ration field, then the density of magnons, or spin flips, 
is dilute. In this case, we can safely use the ladder 
approximation^^ to renormalize the interaction ver- 
tex in a controlled manner. In fact, we strictly speak- 
ing analyze the interactions for fields above the satu- 
ration field, where there are no bosons present in the 
ground state, and we consider just two bosons interact- 
ing pairwise above the vacuum. We require the behavior 
in the limit in which the saturation field is approached, 
i.e. in which the energy of the two interacting bosons ap- 
proaches zero. This limit should be familiar from ultra- 
cold atomic systems, in which the complicated interac- 
tions between atoms can be replaced by one or a few 
scattering lengths, which represent the effective interac- 
tions in the dilute limit. Here we obtain the effective in- 
teractions from the Bethe-Salpeter (BS) equation, which 
reads 



(44) 



Here T(k,k f ;q) is the irreducible four-point interaction 
vertex taken with all external frequencies equal to zero, 
and Q = 2 (ft, — ft sa t) = —2p. The k' are the incom- 
ing momenta and k + q,k' — q are the outgoing mo- 
menta, as shown in Fig. [lOj From this, one obtains that 
ri =T(Q,Q,0) and Y 2 = T(Q, -Q, 0) + T(Q, -Q, -2Q). 
In Eq. p9| ), we introduced a factor of U into the defini- 
tion of V[q) to enforce the spin- 1/2 constraint, which is 
equivalent to taking the limit U —> oo. This limit in the 
BS language, Eq. ( [44]) , pro vides us with an additional 
constraint which reads^^ 



L 



T(k,k';p) 



p e(k + p) + e(k> - p) + ft 



= 1. 



(45) 



Both Eq. (44) and Eq. fl45b can be applied either in two 



or three dimensions, or for the one dimensional TST; 
in the latter case, the integral over p should be re- 
garded as an integral over p x and a sum over the discrete 
p y = 0, 27r/3, 47r/3. Notice that in two or fewer dimen- 
sions, since e{k) ~ k 2 , V(k) ~ 1 near k = 0, the integral 
is at least logarithmically divergent when Q approaches 
zero. This reflects the fact that weak interactions are 
marginally relevant at the zero density fixed point in 
d = 2, and relevant for d < 2. We use this to our advan- 
tage, since we are interested precisely in this limit: the 
singular parts dominate the vertex function as £1 — >• + , 
and we extract these dominant singular terms analyti- 
cally to obtain the asymptotic behavior. For d > 2, the 
integrals become non-singular, and one can directly take 
the H = limit. 



12 



2. Calculation of Y\ and Y 2 in 2d 

We first give a brief summary of our calculations for 
the 2d case. The dispersion minima occurs at k, k 7 = 



±Q 2 d = ±(Q2d,Q2d/2), where Q 2 d is given in Eq. p3 ) 
To solve the BS equation, we use the following ansatz: 



r(fc, fc'; q\ Q) = Aq + A\ cos q x + A 2 sin ^ + As cos g y + A4 sin q y + A5 cos(q y — q x 

I 



A 6 sin(<^ - ) 



(46) 



where are coefficients dependent on fc, fc', J, J 7 and f2. 
With Eqs. (|44||45) |46j , one can solve a set of linear equa- 
tions for the coefficients which gives an explicit form 
of T(q) for a given set of /c, k' ', J, J' and Q. Details of the 
2d case are given in Appendix |B 1| From the solution, 
we simply obtain 



ri >r 2 



for < R < 1, 



which implies that for all range of anisotropics, < R < 
1, the ground state near saturation field is always an in- 
commensurate planar (or fan) state. 

To see how the incommensurate planar state dominates 
over the cone state in the weakly coupled chains region, 
we expand the expression of T's in the leading order of 
both 1/lnfi and j = J' / J 

r,/J= l-^j + lf + oif)]^ 

+[-8jTr\n(4j) + a + 0(f)\ 



'(In ft) 2 



-8j7rln(4j) + 0(i 3 )]- 



a 



'(In ft) 2 
8j7r(24-161n2-37rln2) 



16 + 37T 



> 0. 



(48) 



Since the extra factor a is always larger than zero, the 
ground state always prefers the fan state in the limit of 
decoupled chains. 

One can analytically check this result in the same limit, 
J' <C J. We discuss this extension in Appendix B 3 



3. Calculation of Y\ and Y2 in the TST 

We now present a brief overview of our calculations on 
the TST. We consider an infinitely long system, where 
q x is continuous and q y = 0, 27r/3, 47r/3 is discretized by 
periodic boundary conditions. The dispersion minima oc- 
cur at k,k 7 = ±Qi d = ±(Qi d ,27r/3), given in Eq. ([34]). 
We are now in a position to solve the BS equation, where 
we follow similar procedures as the two-dimensional case. 
We use the same ansatz, Eq. (46), to solve for the coef- 
ficients A*. From these coefficients, we can obtain the 



explicit forms of r(g), for which we provide details in 
Appendix |B 2| Our results are as follows 



ri > r 2 , 
ri < r 2 , 



for < R < 0.48, 
for 0.48 < R < 1. 



(49) 



This tells us that for R < R r = 0.48, the incommensu- 



(47) rate (fan) state is favored, while for R > R C1 the cone 



(umbrella) state is favored. This result is in agreement 
with the analytical result, in Appendix |B 3[ where it was 
shown that spins order into a cone state in the decoupled 
chains limit. 



E. Commensurate-Incommensurate Transitions 
(CIT) 

In the previous subsection, we found that near satura- 
tion, the ground state of the two-dimensional model for 
all R and of the TST for R > 0.48 is coplanar, with mod- 
ulation of the z-component of the spin at wavevector 2Q. 



As mentioned in Section IV C this implies spontaneous 
breaking of the discrete translational symmetry, which 
is sensitive to commensurability effects via the terms in 
Eq. ( plj ). In particular, we expect that the wavevector 
Q will lock to commensurate values, where 2Qn is a re- 
ciprocal lattice vector, over a finite range of field and 
anisotropy, R. We now turn to a description of these 
commensurate-incommensurate transitions (CITs), both 
in the 2d case and for the TST. 

To study the C ITs, we must now consider the full ac- 
tion, Eqs. ( 36|41 ), for h < /i sa t, i.e. for \i > 0, where 
the bosons are at non-zero density. In two-dimensions, 
we can regard them as condensed, while in the TST, 
true condensation is impossible but the system can be 
viewed as a quasi-condensate or a Luttinger liquid. In 
either case, amplitude fluctuations of the ijj a fields are 
small, and we can write the effective action in terms of 
the phases a , where ip a ~ ^ e _26>a in the coplanar/fan 
region. 

Conceptually, the effective action for the phase fields 
is obtained by first following the renormalization of the 
system away from the zero density fixed point, \i = 0, 
where amplitude fluctuations are still important. Once 
the energy scale set by \i is reached, these fluctuations 
are quenched, and it is sufficient to consider only small 
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fluctuations in the amplitudes. To achieve this, we sim- 
ply make the as sumption of small amplitude fluctuations 
in Eqs. ( 36|41 ), but with the bare couplings replaced 
by fully renormalized ones, at the scale p. We believe 
this procedure properly captures the scaling for small /i, 
though it is not quantitatively reliable. 

Because the low energy dispersion of the single magnon 
states is exactly known and described by the quadratic 
terms in Eq. (36), the corresponding couplings are un- 
renormalized. The interactions Y\ and r 2 , however, 
are renormalized by multiple scatterings, which is ex- 
actly what is captured by the BS equation discussed in 
Sec. |IVD| From this analysis, we simply take as our 
renormalized couplings T a (£l — 2p). Note that this 
would be exactly correct if we replaced p by \p\ for the 
case p < 0, but on scaling grounds it should give the 
correct dependence even for p > 0. 

The renormalized interactions can be approximately 
represented for small p as 



r Q (M) 



1 + mu a /((mp) ' 



where 



C(mp) 



(ra/i) 1 / 2 d= 1 

l/|ln(ra/i)| d = 2 



(50) 



(51) 



We can in principle use the renormalize T a (p) for the 
original lattice spin model, which have the same leading 
and first sub-leading terms for small p (up to second 
order in £ <C 1) as in Eq. (50), but with considerably 
more complicated coefficients. Beyond second order in 
C, the lattice T a differ somewhat, and the expression is 
unwieldy. The above form is sufficient for our purposes, 
and is exact for a continuum model. 

Once the T a (p) are known, the analysis 



is 



straightforward 17 . We write ip a = [p + crj 1 / 2 e~ lUa , and 
assume small fluctuations in <j a around the saddle point 
value for 



P 



(ri(/i) + r 2 (/i))' 



(52) 



and u a are constants related to the "bare" values of T a . 



(Here we assume Ti(p) > r 2 (/i)). Eq. (52) properly 
captures, through the dependence of on /i, the non- 
mean-field dependence of the boson density on chemical 
potential. In particular, it yields p ~ p 1 ! 2 in 1+1 dimen- 
sions, consistent with the fact that repulsively interacting 
bosons behave with an effective hard core at low density, 
and consequently have an equation of state similar to free 
fermions. 

Expanding the action to quadratic order in a a and ne- 
glecting irrelevant terms involving derivatives of a a and 
their couplings to higher derivatives of a , we obtain (ne- 
glecting constant terms) 



/ 



d d vdr I ifadrOx + cr 2 <9 T 2 ) 



2m 



|V6»i| 



2 



W) + + a 2 2 ) + rv 1( 7 2 



(53) 



Next, we integrate out the <j a fields, and express the 
resulting action in terms of new linear combinations, 



The result is 



S = Se + Ss, 



where 



with 



and 



2(r 1 ( M ) + r 2 ( /i ))' 



~y^A„ cos[n(9 -q„ -r)]|, 



(54) 
(55) 



So = J d d vdr {^{d T 9f + |(V^) 2 } , (56) 



(58) 



with 



1 



2(ri( M ) 

2w n p n . 



r 2 (M))' 



p 

2m' 



(59) 



Here we have restored the term resulting from S' in 
Eq. (41). Note that the "charge" field 6 describes the 



Goldstone mode of the broken (or quasi-broken in Id) 
U(l) symmetry, and thus remains exactly massless. It 
completely decouples from the field, and can be ne- 
glected in the analysis of the CIT. 

We are no w in a position to analyze the CIT using 
Eqs. ( |58|59 ) and the results of Appendix [A] This is 
strongly dimension dependent, so we treat the cases of 
two dimensions and one dimension separately. 



1. Two dimensions 

In two dimensions, we begin by presuming that one 
of the cosines in Eq. (58) is almost non-oscillating, i.e 
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when one of the q n is close to zero. Generically, this will 
happen for one specific minimal n, when 



Q2d 



Tim 
n 



(60) 



for some specific m,n, with \SQ\ <C 1. The other rapidly 
oscillating cosines can be neglected, and we retain only 
the weakly oscillatory one. Then, in the x, y coordinates, 
the action takes the form given in Eq. ( Al ), with A n = A, 
and q = q n = 2SQ. 

We can now directly apply the re sults of Appendix |A 1| 
Using S = pq = 2p5Q, and Eq. (A6), we obtain that 
the commensurate state is stable for \SQ\ < SQ C , which 
defines the location SQ C of the CIT as 



SQc 



(61) 



where we used Eq. (59) for d = 2, and, of course, we 
assume p > 0. Here 



T(/i) = 



1 



r 1 ( / i) + r 2 ( / i) 



2\ln(mp)\ 
m 



for p < 1, (62) 



is a weak logarithmic function of p. 

For the commensurate state centered around R = 
(J 7 = J), we have n = 3, and the phase boundary for the 
C-IC transition is linear in p, up to logarithmic correc- 
tions. However, as n increases, the widths of the com- 
mensurate phases decrease. 



2. One dimension 

In the TST, to derive the Id theory we must sum over 
discrete y. This restricts the A n terms in Eq. (58) to n 



which are multiples of 3, so that the y component of q n 
(= 2nQ y ) is a multiple of 2tt. 

Following the discussion for two dimensions, we again 
consider wavevectors 



_ Tim 

Qid = h oQ, 

n 



(63) 



with appropriate m, n such that \SQ\ <C 1, and keep only 
the dominant cosine term of order n, which then matches 
the sine-Gordon form in Eq. ( |A1| ) with q — 25Q. Then 
we take over results from Appendix |A2| 

According to that discussion, a commensurate phase is 
stabilized whenever the scaling dimension of the cosine 
term, A n , is less than two. Using the result in Eq. (A12) 
and also, Eq. (59), we obtain 



n 2 fp\ 
V2ti \m) 



1/4 



U\ — U 2 
UiU 2 



(64) 



so that A n <C 1 for p <C 1. This shows that A n < 2, 
and the commensurate phase is indeed realized. Note 
that if we approximate A n = 0, then this becomes the 
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FIG. 11: Entanglement entropy at (a) M/M s = 5/6, = 0.3, 
the incommensurate coplanar phase, and (b) M/M s = 5/6, 
R = 0.66, the cone (or umbrella) phase. We take a system 
size of N x = 120. 



same classical estimate as in the previous section, except 
that T a (p) has a different dependence in one dimension. 
While this is in principle appropriate for very small /i, the 
1/4 exponent in Eq. (64) indicates that A n can be sub- 



stantial nonetheless, so we will proceed with the estimate 
taking A n ^ 0. 

Using S — 2p5Q and the estimate for the critical 5 C in 
Eq. (A16), and applying Eqs. (52) and (59), we find the 



location of the Id CIT as 



SQc - (i 



p 



(65) 



For n = 3 and assuming A n — >• 0, this predicts 5Q C ~ 
/i 1 / 2 , which does not agree with p ~ R scaling of the C-IC 
boundary in the upper left corner of the phase diagram 
in Fig. [3] However the range of p there is not particu- 
larly small, h changes from 4.5 to approximately 3 as R 
changes from to 0.1. This observation calls for a m ore 
careful analysis of behavior predicted by Eqs. ( 64|65 ) for 
p ~ 0(1). We find that numerical coefficients in (64) 
make A n= 3 to vary in the interval 0.5 — 1 for p relevant 
to the C-IC boundary in Fig. [3j resulting in an almost 
linear dependence 5Q C ~ p away from the strict p —> 
limit and in qualitative agreement between our analysis 
here and the numerical data in Fig. [3] 



DMRG results 



In Sec. |IVC| we show that the cone state corresponds 
to a single-Q condensate bosonic field, while the incom- 
mensurate planar state corresponds to double-Q conden- 
sate. This is verified by the central charge measurement, 
where we find c = 2 to describe the coplanar phase as 
shown in Fig. [TT^ i, as opposed to c = 1 for the cone in 
Fig.|TIJ). 
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(2) M/ M a = 5/6, R = 0.30 




10" 



FIG. 12: (Color online) Transverse spin-spin correlation func- 
tion for M/M s = 5/6, N x = 120 at (a)R = 0.30 in the incom- 
mensurate coplanar state and (b)R = 0.66 in the cone state. 
Our DMRG data points are plotted in (black) circles, while 
the theoretical fit, Eq. (66), is shown as a solid (blue) line. 
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The transverse spin-spin correlation function for the 
cone state can be written as 



s\\ / J{e{v)-e{v')) 



(S+S~,) ~ ^ cos(Q.(r-r'))(e 



—2 



^ cos(Q-(r-r / ))C 7? (x,x / ) (66) 

With C v (x,x f ) given in Eq. We fit the DMRG 

results to this formula in Fig. |12]o. The transverse corre- 
lation shows a clear sinusoidal pattern with incommen- 
surate wavevector Q = (1.107T, 2ty/S) and r] = 0.37 at 
M/M s = 5/6, R = 0.66. Fig. [l2]o shows an excellent fit 
which yields the exponent 77 = 0.37. 

The whole procedure is repeated for the incommensu- 
rate planar state, 

(S+S;,) ~ # 2 (cos(Q-r + ^>))cos(Q-r / + ^> / ))) 

9 



■ cos(Q • (r - r / ))C r?+r ^(x, x') 

The exponent 77 and 7) come from averaging the 9 
and fields, respectively. The fitting estimates Q = 
(1.267r,27r/3) and 77 + 77 = 0-54 at M/M s = 5/6, R = 0.3, 
shown in Fig. [12| l 

Next we consider the vector chirality (VC), which is 
defined as V Xi y = z • (S XiV x S x +i^ y ) in Eq. ( |43| ). As dis- 
cussed in Sec. IV C[ since the cone state favors XY order, 
the VC should be a nonzero and constant value. Indeed, 
as shown in Fig.[l3j the VC correlation function does not 
decay with distance in the cone state, i.e., R = 0.66 and 
0.80, and the finite-size scaling (Fig. 13 'b)) shows that 



the corresponding VC order parameter remains finite in 
the thermodynamic limit. Instead, for planar states, the 



FIG. 13: (Color online) Spin vector chirality (VC) correlation 
function, as defined in Eq. (43), at M/M s = 5/6 with system 
size N x = 120. In (a), we show R = 0.4 (purple diamond) 
where the system orders into an incommensurate coplanar 
phase. Furthermore, R = 0.66, 0.80, where the system is in 
the cone phase, is shown on the same plot. We can see that the 
VC approaches a constant in the cone phase while decaying 
in the coplanar phase. In (b), we show the finite-size scaling 
of the VC order parameter in the cone phase. 



spins are confined to one plane, so the VC correlation 
decays exponentially (see R = 0.4 data in Fig. 13). 



V. WEAKLY COUPLED CHAINS 

A. Bosonization of a Heisenberg chain 

In this section, we give a brief overview of apply- 
ing Abelian bosonization to a single spin- 1/2 Heisenberg 
chain in a magnetic field. The Hamiltonian of interest is 
as follows 



H ch = jY, S(x) -S(x + l)-hJ2 S z {x), 



(68) 



x=l 



x=l 



where the magnetic field is chosen along the z-direction, 
and the lattice spacing has been set to 1. Here, the mag- 
netization, M = j-S z (x), is conserved, and hence, 
the magnetic field, /i, can be treated as a chemical po- 
tential to relate the properties at ft ^ to those at 
ft = 0. For any magnetizations less than saturation, i.e. 
M < M sat = 1/2, the low energy theory can be described 
by a canonical set of a massless scalar field, 9, and its dual 
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field 



Hq = J dx V -{{d x cj>) 2 + {d x 9) 2 ). 



(69) 



These two fields satisfy the familiar commutation rela- 
tions 



[0(x),<j>(x')] = -iG(x-x f ) 



(70) 



where 6 is the Heaviside step function. The spin veloc- 
ity, in Eq. (69), is a function of the magnetization, M. 
When M = 0, v/J = n/2, and the 517(2) symmetry is 
restored. For the case when M > 0, v decreases continu- 
ously and is numerically determined by the Bethe ansatz 
integral equations (see Fig. 9 of Ref . [25]) . 

At a fixed magnetization, both the longitudinal (along 
the field direction) and transverse (perpendicular to the 
field axis) spin fluctuations have gapless excitations. The 
longitudinal modes occur at commensurate wave vector 
k x = and incommensurate ones k x = 7r±2J, where 8 = 
7rM, while the transverse modes are at commensurate 
wave vector k x = tt and incommensurate vectors k x = 
±26. Then, one can expand the spin operator around 
these low energy gapless modes, i.e. 



M + S z Q {x) + e^- 2S >S:_ 25 {x) 



S z (x) 

S+(x) = e- m *S± 2S (x)+e i2S *S+(x) 



-25\ 

+(-l) x S+(x), 



(71) 



where Sq, S^ ±2S (x), S± 26 (x) and S n are operators whose 
scaling dimensions depend on M. One can rewrite these 
operators in terms of the bosonic fields, <ft and 0, 



S z {x) = p~ x d x <t>, 



s. 



7T-25 



(X) 



si{x) 



-27Ti0//3 



± l -A 2 e^ e e ±i2 ^l^ 



(72) 



Here, the parameter (3 = 2ir1Z is related to the compatifi- 
cation radius 1Z and can be calculated by solving the in- 
tegral equations, which can be found in Refs. 26-28. The 
compactification radius takes on a simple form, 2ttT1 2 — 1 
at zero magnetization, and approaches 2tt1Z 2 = 1/2 as 
M M sat = 1/2. The constants, A 1 , A 2 and A 3 , 
are determined numerically 11 . Furthermore, at M = 0, 
the scaling dimension of Sq and S± 2§ (x) is 1, and these 
operators can be written in its SU(2) symmetric form 
M = J# + J l • The scaling dimension of S* ±25 (x) and S n , 
however, is 1/2 at zero magnetization and is related to 
the staggered Neel order, N, and dimerization e. Further 
details for the M = case are provided in Appendix |C 1| 
Now, in order to compare our DMRG results to this 
analysis, we must enforce open boundary conditions (BC) 
along the chain direction to mimic DMRG's BC. This 



can be achieved by introducing two additional "phantom 
sites" at x = and x = L + At these positions, 
we enforce boundary conditions on the bosonic field, 0, 
where <j)(x = 0) = and <j){x = L + 1) =0. The sum in 
Eq. (68) now runs from site index to L, and we effec- 
tively obtain a periodicity of L + 1 using these phantom 
sites. We can now substitute Eq. (72) into Eq. (71), and 



enforce the open boundary conditions. The spin opera- 
tors can now be written as (for brevity, we suppress chain 
index y) 



S z (x) = M 
S+(x) = e^^[A 3 (-l) 



^- AlSin( 7^ ) - (7r " 2 ^ ) ' 



2tt 

+A 2 sin(— (t>{x) +26x)], 



(73) 



where M = ML/(L + 1) and 5 = ttM. The bosonic field, 
0, can also be expanded in terms of its lattice modes as 



d>(x) = J2 



sin(q n x) 



(a n + a+), 



(74) 



where q n = 7ra/(L + l). Here, a n and a+ are the annihila- 
tion and creation operators and satisfy the commutation 
relation [a n ,a+,] = 5 n>n /. 



B. Triangular spin tube 

We now extend our previous discussion to study the 
behavior of the TST, described by Eq. ([!]), in the limit of 
weak coupling, J' <C J. Using the low energy expansions 
of the spin operators in Eq. (71), we can express the 
low energy Hamiltonian as H = Hq + -Hi, where Hq 
is described by a sum over the free bosonic modes in 
Eq. (69) on each chain. Here, Hi describes interchain 



interactions and is as follows 



3 ^ 

#i = J 'E / dx{2M 2 + 2S;. fi S; +1 . fi 
v=K=v, 



(75) 



x=0 

Y (i _ e 2icrS )S z i9 ,S* 

\ + e 2iaS\ 



-2a5 



E 

<7=± 



where again, M ML/(L + 1). 

The first term, 2M 2 , with scaling dimension 0, is the 
most relevant, but is trivially a constant. The second 
term is marginal with scaling dimension 2, and renor- 
malizes the Luttinger parameters and the velocities of 
the bosonic fields, 0,0, in Eq. (69). The third term is 
relevant at M = with scaling dimension 1, and be- 
comes marginal as magnetization increases, approaching 
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a scaling dimension 2 as M — >• M sat . This term is re- 
sponsible for the SDW phase that arises when relevant. 
The fourth term, which involves a derivative, is marginal 
at M = with scaling dimension 2 and becomes increas- 
ingly relevant with increasing magnetization, saturating 
to a scaling dimension of 3/2 as M — »> M sat . This is a 
"twist" term that favors the cone or XY phase that or- 
ders perpendicular to the magnetic field. The last term 
is always irrelevant, with scaling dimension > 2 and can 
be neglected in the analysis of this theory. 

Apart from the trivial constant term, the SDW and the 
"twist" terms are the most relevant ones and have com- 
peting scaling dimensions as magnetization varies from 
to saturation. With the exception of some subtleties 
that arise from the TST boundaries (we discuss this 
in proceeding subsections), standard scaling arguments 
can be made about these two operators. For small M, 
the SDW term dominates, and the system orders into a 
collinear SDW in which the ordering momentum, it — 28 
scales linearly with magnetization. The twist interac- 
tion dominates over the SDW at a larger magnetization, 
and the system orders into a cone-like state. Since there 
is no spontaneous breaking of continuous symmetry in 
one dimension, the SDW and cone order are not really 
ordered states, but are Luttlnger liquids with one gap- 
less mode. This competition between cone and SDW 
phase was discussed for 2d triangular lattice in Ref. 0J 
where critical magnetization, M CYltl at which the quan- 
tum phase transition from the SDW to the cone phase 
takes place, was evaluated. The TST has the same criti- 
cal M cr it = 0.64M sat as the 2d case, except that the cone 
state obtained in this quasi- Id regime is smoothly con- 
nected to the cone phase obtained in the high field region 
in Sec. El 



Eq. ( 75 ) is not complete as it does not account for sev- 



eral less-obvious relevant terms which are allowed by the 
lattice symmetry of the problem. This will be consid- 
ered in more detail later. Within the SDW phase, it is 
possible to lock the SDW momentum to a commensurate 
value by accounting for high-order umklapp processes. 
The first of these leads to a commensurate SDW, which 
is in fact identical to the 1/3 plateau with the "up up 
down" structure. This is discussed extensively later in 
Sec. Ell 



Other more relevant intra-chain interaction terms may 
appear due to fluctuations that are not accounted for in 
the naive bosonization in Eq. (75). We will discuss these 
effects in Appendix |C 2| 



lows 



H sd W = ^ / dX~ [(d X (j)y) 2 + (d X 0y) 2 ] + -^d X (t)yd X (j)y +1 

y=i J Z P 



7 sdw cos[ 



,2* 



0y+\ 



■26, 



(76) 



where the bare SDW coupling is given by 7 s( j v 
J'A\ sm{6) > 0. 



1. Scaling considerations 

Renormalization group arguments give considerable in- 
sight into the physics of Eq. (76). All but the last term 
in H s d w are scale invariant, and can be considered a fixed 
point Hamiltonian. The remaining SDW term, propor- 
tional to 7 s dw, is not, and renormalizes under the scale 
transformation x —> bx, according to the usual linearized 
relation 



7sdw0) 



Tsdw? 



(77) 



where b > 1 is an arbitrary scale factor. As discussed 
in the previous subsection, A s d w < 2, so that the SDW 
interact is relevant, and grows in strength under rescal- 
ing. Eq. (77) is valid for small dimensionless 7 s dw(fr)> and 



therefore the weak coupling regime is limited by the con- 
dition 7 s dw(fr) < v. This defines an "SDW correlation 
length" £ sdw such that %dw(b) = v- 



Csdw 



(f/7sdw) 1/(2 - 



v) 



(78) 



In the weakly coupled chain regime, 7 s d w is small and so 
£ s dw is large. On scales large compared to this correla- 
tion length, we expect that the bosonic modes appearing 
inside the SDW term become "pinned" to values which 
minimize this interaction. This pinning corresponds to 
the creation of well-established SDW order. 

Due to the divergence of £ s dw, however, the establish- 
ment of SDW order can be prevented by finite size effects, 
even for reasonably large systems accessible by DMRG. 
For a finite system of length L, we must compare the 
SDW correlation length to L, and it is expected that 
physical quantities will be functions of the dimensionless 
ratio S sdw = £ sdw /L. For S sdw <C 1, SDW-like behavior 
is expected, but when S sdw > 1, there may be a non- 
trivial crossover. This occurs particularly in the case of 
the TST, for which an analysis, detailed below, shows 
that the crossover is discontinuous. 



SDW 



2. L = oo 



In the region of low to intermediate magnetization 
and small J' , we can neglect all terms in Hi except the 
marginal one and the SDW interaction. Using bosoniza- 
tion, Eq. (73), the Hamiltonian can be re- written as fol- 



For an infinitely long system, 5 = 0, we can under- 
stand the nature of the SDW state by simply minimizing 
the 7 s dw term in Eq. ( 76 ). When the width is also infinite, 
i.e. in two dimensions, one can simultaneously minimize 
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each cosine term (for each y) independently. This occurs 
by taking 



2tt 

T 



25 



-y, 



(79) 



d=2 



where tp is an arbitrary constant (x- and ^/-independent) 
phase. Allowing for small gradients of (p, which might 
be present due to fluctuations or perturbations, the spin 
operator can then be represented as 



S*(x)\ 



d=2 



M 



2tt 



A 2 sin [p( x )-7^l(2x-y)], 

which indeed is the classic form for a spin density wave 
with wavevector n ~^ s (—2, 1). This corresponds to an 
ideal two dimensional SDW state, and cp gives the "slid- 
ing" or "phason" mode of the SDW. For generic irra- 
tional 5/tt, remains a gapless pseudo-Goldstone mode 
associated with translational symmetry breaking. In two 
dimensions, the zero point fluctuations of this mode do 
not, however, destroy long-range SDW order. 

Now consider the case of the TST ladder, where y = 
1,2,3 and periodic boundary conditions are applied. In 
this case it is generically impossible to simultaneously 
minimize each cosine term separately. Instead, the min- 
imum occurs when 



2tt 

T 



2tt 



-y, 



(81) 



L=oo,TST 



where again cp is an arbitrary constant, reflecting the 
invariance of Eq. (76) under uniform translations of all 
the (j) y . Again, one can express the spin operator here 
using this form 



L=oo,TST 



M 



2tt 

-A 2 sin [(p(x) 



(82) 

(7r-2£)x+f y]. 



In contrast with Eq. (79), the minimum configuration in 
the TST, Eq. (81) is independent of 5, manifesting in 



Eq. (82) as a difference dependence on y from Eq. (80). 
The difference is due to the frustration of the intrinsic 2d 
SDW order by periodic boundary conditions, which tend 
to lock the SDW order to a commensurate form in the 
y direction. Interestingly, the two results coincide when 
5 = 7r/6, which corresponds to the case M = M sat /3. 
At this point, the periodicity of the TST and the SDW 
order are compatible. 

As in the 2d case, at the level of Eq. (76) applied to 
the TST, the uniform translation mode cp remains gap- 
less. Unlike the 2d case, however, in one dimension, the 
zero point fluctuations of this mode are sufficient to dis- 
rupt long range SDW order, which instead manifests as 
power law correlations. Nevertheless, the short distance 
physics is still that of an SDW, and moreover the Id 
fluctuations are easily accounted for theoretically. This 
is accomplished simply by treating cp as a free massless 



boson, as we discuss below in Sec. V C 3 



3. Finite length L < oo 



As we have discussed in Sec. |V A[ for a finite length 
chain, we must impose the boundary conditions (j) y (x = 
0) = 4> y (x = L) = 0. These conditions are incompatible 
with the values, in Eq. (81), which minimize the SDW 



term in the infinitely long case. This means that end ef- 
fects strongly affect, and tend to suppress SDW ordering. 
What do we expect? For short systems, where S ^> 1, 
the end effects will dominate, and the effects of the SDW 
interaction become negligible. In other words, all com- 
ponents (j) y will be largely not affected by the SDW term, 
and the system should behave similarly to three decou- 
pled chains of finite length. For long systems, S <C 1, the 
SDW pinning should be effective far from the boundaries, 
and only the pseudo-Goldstone mode $o will behave like 
a massless field (pinned at the boundaries). 

Let us now address the crossover. It is convenient to 
first make a change of basis^ from the <pi , fa •> 4>3 to new 
fields $o, $2: 

(l/y/3 l/y/2 l/y/6 
n I = 1/^3 -2/V6 
\l/V3 -l/y/2 l/y/E 

The dual fields 6 y transform similarly. Note that the cen- 
ter of mass field is just proportional to the SDW phase 
introduced earlier: <£ = ^^ ( P- The boundary condi- 
tions (j) y = at the ends translate to ^ = at the ends. 









• (83) 







The SDW Hamiltonian now reads H S( ± W 
where the harmonic part 
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sdw 



■H. 



(i) 

sdw' 
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sdw 
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n=l J 



dx 



2k, 



(d x e n f 



(84) 



is expressed in terms of renormalized stiffnesses ft 2 = 
1 + 4J'/(/3 2 v) and = 1 - 2J' /(j3 2 v) and velocities 
v n = v/n n . Its interacting part (the analog of the second 
line in Eq. ([76]) written in the new basis) reads 



^sdw — 



/• 



" sdw / dx 2 cos 
cos[-^-V2^>i 



7T - 25 



$1 



7T 



-25, r 27r 



$ 2 



(85) 



Note that the center-of-mass mode &o ex (p does not enter 
in Eq. (85). Thus it behaves as a free massless boson, 
independent of the strength of the SDW coupling. The 
distinction between S and S in the SDW Hamiltonian is 
not important when analyzing the crossover, and will be 
dropped in this subsection from now on. 

To analyze the crossover, we first carry out the renor- 
malization group procedure by integrating out fluctua- 
tions of the fields due to modes with wavelength less 
than the system size L. In doing so, we replace 7 s d w 
by its renormalized value at this scale, 



Tsdw 



— )> 



7sdw(£) = L Asdw 7sdv 



(86) 
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Note that we have done the coarse-graining step of the 
RG of integrating out modes, but we have not rescaled 
any fields or coordinates, so as to keep the original units 
unchanged for clarity. Under this coarse-graining trans- 
formation, the quadratic terms in the Hamiltonian re- 
main unmodified. 

In this renormalized Hamiltonian, it is appropriate to 
carry out a classical saddle point approximation for <3>i 
and $2, which are the fields pinned by the SDW coupling. 



The SDW potential in Eq. (85) is minimized by $2 = 0, 
which is compatible with the boundary condition, and 
so, we can impose this condition. Then only $1 enters 
the saddle point condition in a non-trivial way. For sim- 
plicity we specialize to the case £ = 7r/6, orM = M sat /3. 
Then we may define \I/ = ^f^i + x> f° r which the sad- 
dle point Hamiltonian, neglecting the decoupled $0 term 
becomes 



class 



y 4L ^{K(S^) 2 - 7 sdwW(cos[2^+2cos[^])}, 

(87) 

with K = /3 2 £i/47r 2 fti. 

The 7 s dw term is clearly minimized by \I/ = 0, while the 
open boundaries require ^(0) = ^(L) = 27r/3, causing 
the strong suppression of SDW order by the ends. There 
can be a non-trivial configuration, Sfr(x), which minimizes 
the functional H c \ ass . To bring out the crossover physics, 
we transform to dimensionless coordinates, letting 



x ^K/% dw (L)z, 



which gives 



with 



L 

£sdw 



V K lsdw(L) 
(L/0 1-A sdw 

(^/7sdw) 1/(2 " 



/2 ^A sdw /2-l 
"sdw ' 

v) 



(90) 
(91) 
(92) 



Note that Eq. ( |92[ ) agrees, at the level of scaling, with 
Eq. (78) obtained earlier from general arguments. For 
the purpose of minimization, the overall prefactor eo is 
irrelevant, so it is clear already from Eq. (89) that the 



properties are a function of the scaling variable £ s dw only, 
as expected. 

We are now prepared for the saddle point approxima- 
tion, which consists in minimizing Eq. (89). Starting 



from the Euler-Lagrange equation, which has the usual 
"energy" integral of motion, one obtains 



dz 



= C-(2cos[^]+cos[2^]), 



(93) 



where the integration constant ("energy") C is fixed by 
the condition d^(z = L/2)/dz = as C = (2cos[^i /2 ] + 
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FIG. 14: (Color online) S z profile for R = 0.7 at (a) M/M s = 
1/3 plateau and (b)M/M s = 1/2 in the SDW state for the 
non- frustrated chain (see text). We show the DMRG results 
in (black) circles and the theoretical prediction, in the con- 
tinuum limit in Eq. (95) in (red) solid line/square. 



cos[2^ 1/2 ]), where we denote \P 1/2 = V(z = L/2). As 
a result the mid-ladder value of \I/ is implicitly given by 
the following integral 



L /-27T/3 

H class = e ^{(S^) 2 -(cos[2^]+2cos[^]) , (89) / 

JO L J J *l/2 



dtp 



* 1/2 y/2 cos[^i /2 ] + cos[2^ 1/2 ] - 2 cos[. 



• cos[2<^] 
(94) 



The full crossover (in this saddle point approximation) 
is obtained from Eq. (94). First, we observe that in the 



limit ^1/2 — > 0, the above integral diverges logarithmi- 
cally, implying that indeed, ^(L/2) = in the infinite- 
size limit. The short system size limit is less obvious. 
For small L, we must choose \&i/ 2 to minimize the in- 
tegral. However, if we make the obvious choice to let 
^1/2 = 27r/3 — e, with e — » + , one finds that the integral 
in fact does not vanish but approaches the constant value 
ir/y/6. In fact, the integral as a function of ^1/2 has 
a non-monotonic dependence, and the minimum value 
of the integral is « 1.1436 < tt/VQ = 1.2826, which is 
achieved for tf 1/2 « 1.3178 < 2tt/3 = 2.0944. Regard- 
less, the lower bound on the integral implies that there 
is a minimum dimensionless length, L m i n > 2.28, such 
that for L < L m i n , the minimum action solution is sim- 
ply ^ij 2 = 27r/3, i.e. V(z) = 1/2 for allz. For such short 
systems, the boundary conditions completely disrupt the 
SDW order, and the system behaves as though it were 
just decoupled chains. The transition from L < L m - m 
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FIG. 15: S z profile of DMRG result for SDW state at 
M/M s = 1/2, R = 0.7 for frustrated chains, (a) y = 1 and 
(b) y — 3. We can see from these plots that the transla- 
tional symmetry is spontaneously broken, and that the SDW 
is strongly affected by boundaries. 
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to L > L m i n is evidently discontinuous, since ^1/2 must 
jump from a value ^1/2 < 1.3178 at L = L m i n + e to 
^1/2 = 27r/3 for shorter systems. To precisely determine 
the value of Z/ m in requires a comparison of the action 
of the non-trivial and trivial solutions to see where they 
cross. 

What are the consequences of this transition? In nu- 
merics, the transition can be probed by varying L or 
varying J' / J at fixed L. In either case, on crossing the 
transition, one expects a sharp change from SDW-like be- 
havior for L > Z/ m in to decoupled chain-like behavior for 
L < L m - m . In the SDW-like regime, the two modes 3>i, $2 
may be considered to have developed a gap, and conse- 
quently, the entanglement entropy of a bipartite cut of 
the sample is reduced compared to the decoupled chain- 
like regime. Specifically, in the SDW-like regime a log- 
arithmic growth with L is expected and consistent with 
central charge c = 1, while in the decoupled chain regime, 
the behavior should be closer to c = 3. At the transition, 
a sharp drop with increasing L of the entanglement en- 
tropy is expected. More detailed predictions can be made 
for the spin density profile, (Sy(x)). We make such a 
comparison in the following subsection. 



D. DMRG results for SDW 

A number of measurements in the DMRG give evi- 
dence of the SDW state. As discussed in the previous 
subsection, the SDW regime of long TSTs can be de- 
scribed by pinning the fields $1 = $2 = 0, and al- 
lowing for gapless fluctuations of the free massless bo- 
son field In the semiclassical approximation dis- 
cussed in Sec. |VC 3[ one can do somewhat better by 
using the $0 fluctuations and replacing <£>2 —> and 

$i(x) ffipfrpE) - if), with ifr(x) given by the so- 
lution of Eq. ([93]) . In this way, one obtains from Eq. ( 73 ) 



(Sfc)) = M 



y 



2tt 



(95) 



FIG. 16: Entanglement entropy for SDW phase for R — 0.5 at 
(a)M/M s = 1/6 below the 1/3 plateau and (b)M/M s = 1/2 
above the 1/3 plateau. Due to large finite-size effects of this 
measurement, we chose to run our simulations on a larger 
system size, N x = 180, in (a), to compare to a smaller size, 
N x = 120 in (b). 



Here the quantity 



X = [^±Hsin(^i 



(96) 



arises from the quantum average over the free boson field 
$0, which is evaluated along the lines of Ref. 11 1 with the 
result that the exponent 



Vsdw 



7T^o 

3/3 2 



kq 1 
6 2ttK 2 ' 



(97) 



For M — M sat /3 and small J', we estimate ^0 ~ 1 and 
2ttK 2 « 1 - 1/(21^6^8/(7^)]) = 0.72 (see Appendix 
A of Ref. 4), which leads to 7^ s dw ~ 0.23, so the spin 
density profile decays quite slowly with distance from the 
boundary in the SDW regime. Note that the y = 2 chain 
does not depend on so one can directly compare the 
numerically obtained magnetization profile for the 'non- 
frustrated' chain with Eq. (95), see Fig. 



14 



2tt 



below. 

One may wonder about the selection of the y = 2 chain. 
For the geometry of our simulations, the model has full 
translational symmetry, y y + 1 in the y direction. 
This symmetry is broken by our combined choice of sad- 
dle point \£ = $0 = in the bulk and the boundary 
condition $0 = at the edges. Examination of the inter- 
action term in Eq. ( 85 ) shows that there are apparently 
two other minimum solutions, ^ = tt and $2 = =L/3/>/6- 
In the infinite system, these are equivalent to the one we 
have chosen, insofar as they give identical results for all 
operators if we make a suitable translation of <£>o. How- 
ever, the choice of boundary condition for $0 prevents 
this translation and results in a broken symmetry state. 
By a different choice of the otherwise equivalent saddle 



sin[(2 -»)(*(*)-— )-(tt-2<5)4 

points, we can obtain formulae analogous to Eq. (|95|) but 
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FIG. 17: (Color online) The transverse spin-spin correlation 
function on a log- linear scale for R = 0.7 and at magnetiza- 
tions (a)M/M s = 1/6 and (b)M/M s = 1/2, both in the SDW 
state for system size N x = 120. Data points are shown as 
(black) circles while the (red) line is a fit to a pure exponen- 
tial function. 



FIG. 18: (Color online) The "octupolar" correlation function, 
((Ul =1 S+(x))(ILl =1 S~(0))), shown on a log-log scale. We 
show for R = 0.7 at magnetizations (a) M/M s = 1/6 and 
(b)M/M s = 1/2 in the SDW state. Our DMRG data points 
are p lotted in (black) circles, while the theoretical fit to Eq. 
(|100| is shown as (blue) line. 



with the y = 1 or y = 3 chains independent of In 
principle, for a finite system even the discrete transla- 
tional symmetry should be unbroken, but the restoration 
of this symmetry is probably only at extremely low en- 
ergies at which tunneling occurs between these minima, 
and indeed we find the symmetry to be spontaneously 
broken in our DMRG simulations. 

In the decoupled regime, L <C -^min? if is more appro- 
priate to just calculate the spin expectation value using 
the free theory, Eq. ([84]), for all three fields $ , $2- 



Then we obtain, instead of Eq. (95), the result that 



(S*(x))=M+^-sm[(Tr-2~5)x}, (98) 
where the "decoupled chains" exponent is 



7r(ft ± fti ± ft 2 ) 

W 2 " 



(99) 



In the same small J' approximation, this gives rjdc ~ 
3?7sdw> so that 77dc ~ 0.610. Note that there is a much 
more rapid decay of the spin density profile from the 
boundary in this regime. 



We compare the spin density profile in Eq. (95) with 



our DMRG data and find reasonable agreement. Fig. [14] 
shows a comparison of numerical data with magnetiza- 
tion profile of the non-frustrated chain, i.e. the y = 2 
result of Eq. (95), while Fig. 15 shows that of frustrated 



chains, y = 1,3. 

We can also measure in DMRG the central charge via 
entanglement entropy, which yields c — 1 for the SDW 
phase as opposed to c = 3 for decoupled chains. This 
is shown in Fig. [l6j where the plots show that at mag- 
netizations M/M s = 1/6,1/2 for R = 0.5, the central 



charges obtained from numerics are c = 0.9, 0.95, respec- 
tively. These values are very close to the predicted c = 1, 
which gives evidence for the SDW. 

Another measurement we can perform is the transverse 
spin-spin correlation function, which should decay expo- 
nentially to support the SDW state. We observe exactly 
this behavior from our simulations, as shown in Fig. [l7[ 
Finally, power-law behavior is expected for the "octupo- 
lar" correlation function^, 

((Ill =1 S+(x))(Ul =1 S-(x'))) ~ C m {x,x'). (100) 

The operator Uy =1 S+ (x) may be though of as inserting 
a soliton - an extra period - into the SDW. This corre- 
lation function decays in the thermodynamic limit with 
the power-law exponent 



3f3 2 

27Tftn 



(101) 



We indeed observe such power 
DMRG, as shown in Fig. [l8 
M/M s = 1/2, R = 0.7) gives 773 = 



law behavior in the 
Fitting this data (for 
3.09±0.2, while the S z 



profile in Fig. [14] for the same parameters is fit to ?7 s dw = 
0.2 ±0.1, yielding the product r^sdw = 0.62 ±0.37. The 
uncertainties for each exponent is crudely estimated by 
tracing out the boundary values when the fitting starts 
to mismatch the DMRG result. The slow decay of the 
S z profile and strong boundary effects as seen in Fig. [14 
induce significant uncertainties in the estimate for ^ s d w , 
so we consider the degree of agreement to the expected 
value r]sr] s dw = 1/2 satisfactory. 
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VI. M = M sat /3 PLATEAU 



Plateau states from SDW 



Magnetization plateaux are observed frequently in 
models of frustrated magnetism, and in a number of 
experiments on such materials. Theoretically, we de- 
fine a magnetization plateau as a ground state of a spin 
system in a magnetic field ft, such that for a range of 
fields, fti < ft < ft2, the magnetization (along the field) 
M(h) = M p is constant. This implies that the magneti- 
zation is a good quantum number, and, since by assump- 
tion the only term in the Hamiltonian coupling to the 
applied field is ftM, that the ground state wavefunction 
itself is independent of the field in this range. More- 
over, since the magnetization M is just the total spin 
S* ot along the field direction, the symmetry under rota- 
tions generated by S^ ot is unbroken. Thus, there can be 
no spin expectation values normal to the field. Further- 
more, no other nearby states must cross the ground state 
(in energy) in this field range, since it remains the ground 
state, and thus, since states with different magnetization 
must have energy depending linearly on the field, there 
must be a spin gap to excitations which carry non-zero 
spin S z relative to the plateau state. 

There are restrictions on such gapped states, follow- 
ing from the Lieb-Schultz-Mattis theorem and related 
arguments 31 . One way to understand them is to map 
the spins to hard-core bosons, where the boson number 
rii = Sf + 1/2. A gapped, insulating ground state of 
bosons in one dimension must have an integer number 
of bosons per unit cell. This implies that the total spin 
XliGu c P er un ^ ce ^ mus t be an integer if the unit 
cell contains an even number of sites, and must instead 
be a half integer if the unit cell contains an odd num- 
ber of sites. Often, such gapped plateau states may be 
considered as ordered states with spins arranged in some 
pattern parallel and antiparallel to the field within a unit 
cell. 

A prominent feature in the phase diagram we obtain 
is a magnetization plateau at one third of the saturation 
magnetization, M = M sat /3. This has been exten sively 
studied in the literature for the isotropic model ! 32 * 33 !, 
R = 0, where it is usually regarded as a result of quan- 
tum "order by disorder". The structure of the plateau 
state in that case is indeed in agreement with a semi- 
classical approach 2 -!, and has a unit cell consisting of two 
up and one down spin, forming a three-sublattice enlarge- 
ment of the primitive triangular lattice unit cell. Based 
on a combination of our DMRG studies and an analytic 
analysis of the quasi-ld limit, J' / J <C 1 (below), we 
show that, in the 2d system, the plateau state persists 
in the full range of anisotropics < R < 1 and forms a 
single phase throughout. For the one-dimensional TST, 
however, we find that the plateau, while present in the 
isotropic regime, terminates before reaching the decou- 
pled chains limit. Both these results can be understood 
from the relation between the plateau state and the SDW 
phase, as will be explained in the next subsection. 



The collinear SDW state shares many of the expected 
elements of the plateau phase. It has an unbroken U(l) 
symmetry, even in the 2d limit, and exponentially de- 
caying transverse correlations in the TST. It has rather 
long-range oscillating correlations of the component of 
the spin parallel to the field, and consequently a markedly 
modulated (Sy(x)} profile in finite systems. The distinc- 
tion between the SDW and the plateau phase is that the 
former is generically incommensurate and gapless. 

Both these differences may be removed due to further 
interactions neglected up to now, which pin the gapless 
phason mode cp at specific discrete values. This has been 
discussed at length already in Ref. 4 for the two dimen- 
sional case. There, it was argued that an infinite sequence 
of plateaux occur at T = within the SDW phase, the 
strongest of these being the 1/3 plateau, and that all 
these plateaux exist at arbitrarily small J' / J . In the two- 
dimensional system, the plateau width (in a magnetic 
field) can be estimated to scale as J(J'/J) 9 / 2 , see Ref. |4j 
Here, we will restrict the discussion to the TST, and find 
that one-dimensional fluctuations suppress most of these 
plateaux, including the 1/3 plateau for sufficiently small 
J'/J. 

The plateau formation is due to additional interactions 
neglected in the sine-Gordon Hamiltonian presented so 
far in Eqs. ( |84|85 ), which involve higher harmonics of the 
phason mode cp. The allowed terms are obtained directly 
from a symmetry analysis. The action of the symmetries 
of the problem on cp may be understood directly from 
the expression for the spin operator in the SDW phase 
of the TST in Eq. (82). Under each symmetry, which 



is a lattice space group operation, cp must be chosen to 
transform appropriately so that Sy(x) is a scalar. This 
dictates the following transformation rules 

1. translation along x, x —> x + 1: cp ^ cp + n — 28. 

2. translation along y, y y + 1: cp ^ cp — 2ir/3. 

3. 2D inversion, x — ^ —x, y —> 2 — y: cp — » — tt/3 — cp. 

In addition, there is a "gauge invariance" arising because 
of the ambiguity of cp due its definition as a phase vari- 
able, which forces the invariance of the Hamiltonian un- 
der local shifts of cp by 2ir. Note that in this section, 
we always consider the infinite L limit, and neglect the 
difference between S and 5. 

Using the local gauge invariance, we seek terms of the 
form 

Hpin = ^ / dxt n sm(mp + a n ), (102) 

n J 

where t n and a n are arbitrary parameters. (In general we 
can also allow a n to be a slowly varying linear function of 
x, which is important for a full analysis of commensurate 
to incommensurate transitions, but we do not require this 
here for the more limited purpose of just identifying the 
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FIG. 19: Entanglement entropy at (a)R = 0, (b)R = 0.2 and 
(c)R — 0.4 in the 1/3 plateau state. We see that the en- 
tanglement entropy approaches a constant for large x , which 
corresponds to a central charge of c — in the ordered state. 



relevant plateau states). Using the translational symme- 
try along y, we immediately obtain the constraint that 
t n = unless n is a multiple of 3, and so we set n = 3k. 
The inversion symmetry then forces a n = (mod 27r), 
so finally, we find 



^pin = ^2 I dxtk sin. Skip, 



(103) 



where we have redefined the tk appropriately. Now it 
remains to apply translation symmetry along x. This 
simply gives the condition that 3k(ir — 25) is an integer 
multiple of 2tt. Writing 5 ttM (7r/2)M/M sat , we 
have 



M 



sat 



3fc - 2p 
3k ' 



(104) 



with k,p integers. This gives a rational family of po- 
tential magnetization plateaux, whose strength decreases 
with increasing k. 

An actual plateau occurs for a given value of magneti- 
zation characterized by integers k,p only if the associated 
term, is relevant 21 , when considered as a perturbation 
to the low energy Hamiltonian of the SDW state, which 
is just the free massless field the ory f or cp. The scaling 
dimension of the operator in Eq. ( 103 ) is easily obtained 
as A 3k = 9k 2 rj sdw = 37r/c 2 ft //# 2 , c.f. Eq. ([97]), and there- 
fore, under RG, we find 



tkip) = t k b 



2-a 3 



t k b 2 



-9/c 2 r7 s , 



(105) 



Here, tk is relevant, and a magnetization plateau appears 
when As/c < 2. Consider the case k = 1, which corre- 
sponds to t he ca se M = M sat /3, and small J' /J. There 
(recall Sec. |VD|) rj sdw w 0.23 so A 3 w 2.07 > 2, and 



FIG. 20: (Color online) Log plot for transverse spin-spin cor- 
relation function at (a)R — and (b)R — 0.2 in the 1/3 
plateau state for symtem size N x = 120. Data points are 
shown as (black) circles while the (red) line is a fit to the 
exponential function. 



with /c, clearly all other potential plateau with larger k 
are absent in the quasi- Id limit. Thus we expect that for 
J' I J <C 1, the SDW state remains stable, and there are 
no magnetization plateaux. 

With increasing J 7 , however, r]^^ decreases, owing to 
its dependence on ^0 in Eq. (97). Including this depen- 



thus t\ is irrelevant. Because A3& increases quadratically 



dence, and using the quasi- Id formula for k (in the text 
following Eq. (84)), we obtain the condition that t\ be- 
comes relevant, i.e. A3 < 2, when J' / J > 0.17. We 
believe that this is still in the domain where the quasi- 
ld approach is valid. The result predicts that the 1/3 
plateau appears only for R < 0.83 in the TST. At fixed 
M = M" sa t/3, the transition from the gapless SDW to 
gapped plateau state at this value of R or J' / J is in the 
Kosterlitz-Thouless universality class, as is well-known 
for the quantum sine-Gordon model. Consequently, the 
gap vanishes exponentially on approaching the transition 
from the more isotropic side, and the ground state energy 
itself shows only an unobservably weak essential singu- 
larity at the transition. We note that other potential 
plateaux with n = 3k > 6 are so strongly suppressed by 
fluctuations that we do not expect any to occur, at least 
in the quasi- Id regime. 

It is interesting to consider the spin structure on the 
plateau. This depends on the sig n of t = t\. For t > 0, 
the sin3(/? pinning term in Eq. ( 103| ) is minimized by 
three values with equal energy, cp = —tt/6 + 27rn/3, with 
n = 0,1,2. For these values, using Eq. ( |82| ), the spin 
density profile takes the form 

(Sfa)) t>0 = M + A X sin [f + f (X - y - n)] . (106) 

This equation describes a three sublattice structure with 
two spins "up", i.e. with (Sy(x)} > M, when x — y — 
n = 0, 1 (mod 3) and one spin "down", when x — y — 
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FIG. 21: £ 2 profile for the 1/3 plateau state at (a)R = 0.0 and 
(b)i? = 0.2. We show data for the three sublattices, showing 
the up- up-down structure, as square (black), circle (red) and 
triangles (blue). There are slight boundary effects near the 
end of the chains, near x — 0, 120; however, the chains are 
well-ordered toward the center, (c) shows the S z profile at 
M/M s = 1/3 as a function of R. We observe that the width 
of the plateau decreases and eventually vanishes near R w 0.8. 



n = 2 (mod 3). This is the semi-classical up- up-down 
state, and has precisely the same qualitative structure as 
predicted semiclassically in the isotropic limit J' = J. 

For the other case, t < 0, the minima occur for cp — 
+7r/6 + 27rn/3, and the spin density profile becomes 

(S z y (x)) t<0 =M-A, sin [f - 2 f(x - y - n)] . (107) 

This describes instead a three sublattice structure with 
two spins nominally "down", with x — y — n = 
0, 2 (mod 3), and the remaining one up. This state does 
not have a natural semiclassical picture, and instead cor- 
responds to the 'quantum' version of the plateau, dis- 
cussed for the two-dimensional lattice in Ref. [4j A car- 
icature of this state is a three site unit cell with two 
sites forming a spin singlet entangled pair, and the third 
(the "up" site) polarized along the field. Our DMRG re- 
sults are consistent with the up-up-down configuration, 
Eq. (|1Q6[), suggesting that t > case is realized. 
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FIG. 22: Finite-size scaling of the boundaries of the 1/3 
plateau with different anisotropics, i.e. (a) R = 0.0, 0.2, 0.4 
and (b) R = 0.6,0.7, 1.0. (c) Width of the 1/3 plateau as a 
function of R. 



We should stress that, apart from the quantitative es- 
timate of fto, nothing in this subsection depends upon the 
quasi-ld approach. The conditions for the existence and 
stability of the plateaux arising out of the SDW state 
are otherwise completely general results based only on 
symmetries of the TST and general arguments. 



B. DMRG results for plateau 

In this section, we discuss how we use DMRG to probe 
into the 1/3 plateau. The first observation of its existence 
is the constant entanglement entropy for the ranges of R 

This shows an 



on the 1/3 plateau, as shown in Fig. 19 
ordered state which corresponds to central charge c = 0, 
in Eq. (12). Furthermore, we can measure the transverse 



spin-spin correlations, which should decay exponentially 
in the 1/3 plateau. We show this measurement in Fig. 20 , 
for R = 0.2 as well as the isotropic case R = 0. 

In Fig. 2lja,b), we plot the S z profile of the spins form- 
ing the three sublattices on the 1/3 plateau. Near x = 
L/2, we see a perfect up- up-down structure, with some 
boundary effects on the edges of the chain. This gives 
definitive evidence of the robustness of the 1/3 plateau 
in these ranges of anisotropies. Moreover, in Fig. 21 ^c), 
we see that the plateau persists up until R w 0.8, at 
which point, the system undergoes a Kosterlitz-Thouless 
transition that destroys the 1/3 plateau. As described in 
the previous subsections, this is a signature of the Id TST 
only: in 2d, the plateau is even more robust, extending 
down to R = 1. This is further discussed in Sec. IVIII Al 

To characterize the properties of the plateau as well 
as its width, we will adopt the following method, which 
takes advantage of the total spin conservation due to the 
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presence of the U(l) symmetry with a magnetic field 
along z-axis. In this case, we can work in a given to- 
tal spin sector S z = J^^f, an d get the corresponding 
ground state energy E(S Z ) 



E(S z ,h)=E(S z )-h-S z . 



(108) 



Then the energy difference between two adjacent spin S z 
sectors is given by 



5E(S Z , h) = E{S Z + l,h)- E(S Z , h). 



(109) 



Generally, at small magnetic field ft, E(S Z + 1, ft) > 
E(S z ,h), so 5E(S z ,h) > 0. However, E(S Z + l,h) < 
E(S z ,h) when ft is large enough, so 5E(S z ,h) < 0. 
Therefore, the boundaries of the plateau can be de- 
termined when E(S Z + l,ft) = E(S z ,h), with the up- 
per boundary h 2 c (S z ) and lower boundary h} c (S z ) of the 
plateau given by 



ft^S 2 ) = £(S* + 1)-£(S*), 
ftJ(S 2 ) = E(S Z )-E(S Z -1). 



(110) 



Finally, the corresponding width of the plateau can also 
be obtained as 



W(S') = h 2 c (S')-hl(S'). 



(Ill) 



In DMRG, the bou ndar ies of the 1/3 plateau can be 
computed using Eq. (110) by fixing the total spin to 
S z = ^jf^- Here, M s = - is the saturation magneti- 
zation, and N is the total number of sites. As shown in 



Figs.|22[a,b), both the upper and lower boundaries of the 
1/3 plateau are determined using different system sizes 
and anisotropics. The corresponding width of the plateau 
is also given in Fig. 22 'c) using Eq. (111). From this, we 



can see that the plateau is very robust and remains finite 
when the anisotropy, R, is small, and decreases with in- 
creasing R. Interestingly, the plateau still remains finite 
even R is very large, i.e., R = 0.7, although the width 
W is very small. In the region 0.7 < R < 1, finite-size 
scaling of the data shows that the width of the plateau 
is zero within the numerical error, for example, at the 
decoupled chain limit R = 1. 



VII. LOW FIELD REGIME 

At zero field, there is already considerable work 
on the spati ally anisotropic Heisenberg model in two 
dimensions^MMl Away from the quasi- Id region, i.e. for 
< R < 0.8, the ground state of the 2d model is unam- 
biguously magnetically ordered, in a coplanar spiral with 
an incommensurate wavevector that varies continuously 
with R. With increasing anisotropy, the ground state 
is less clear, and is quite difficult to resolve numerically, 
owing to the fact that correlations between chains set in 
only at extremely long length scales for small J' / J . A 
controlled renormalization group approach predicts, how- 
ever, that in the limit < J' / J <C 1, the system devel- 
ops a collinear magnetic state instead of the spiral one^. 



Such a collinear state is qualitatively distinguished from 
the spiral one by its pattern of symmetry breaking, which 
leaves a residual U(l) spin rotation symmetry about the 
ordering axis, in contrast to the spiral state which fully 
breaks SU(2) symmetry with no residual continuous in- 
var iance remaining. 

Here we turn to the situation in the one-dimensional 
TST. We argue that in this case the spiral order is con- 
verted by Id quantum fluctuations into a fully gapped 
state with spontaneous staggered dimerization. The ar- 
gument is quite general and is expected to hold for any Id 
system with local non-collinear order and a half-integer 
spin per unit cell. Furthermore, specifically for the TST, 
we show that the tendency to short-range spiral order is 
more robust than in 2d, and unlike in 2d, it prevails over 
collinear order even in the limit of arbitrarily small J' / J . 
Thus staggered dimerization is predicted at zero field for 
all < R < 1 for the TST. See Appendix [CT] for an 
alternative calculation that leads to the same conclusion 
as the one presented below. 

Given the presence of dimerization in zero field, we 
can discuss the behavior in low fields, or more properly 
for small magnetization, in terms of the elementary ex- 
citations of the symmetry broken dimerized state, which 
are domain wall solitons. We obtain in this way differ- 
ent gapless phases at low field, including the SDW state 
discussed previously from the quasi- Id point of view. 



A. Zero field dimerization from spiral order 

In the following, we assume that on short space and 
time scales, the spins establish a similar spiral order to 
that of the 2d system. This notion can be made more 
systematic by considering spin tubes made by wrapping 
the triangular lattice into cylinders with larger circumfer- 
ence. Once the circumference is large enough compared 
to the correlation length of the spiral order, the latter 
should become well-established. It seems reasonable to 
regard this as being the case already for the circumfer- 
ence three TST studied here. This is corroborated also 
by the close correspondence of the phase diagram in the 
weakly anisotropic limit, R <C 1, and the expected semi- 
classical one, as discussed already in Sec. |III| 

With this assumption, the description of the TST 
should be that of a Non-Linear a-Model (NLaM) for the 
spiral order, confined to the finite width cylinder. This 
starting point is similar to the one of Haldane 40 applied 
to unfrustrated spin chains of spin 5, which locally es- 
tablish collinear Neel order. From this formulation, Hal- 
dane established the existence of a featureless gapped 
state for integer 5, while it is known that chains with 
half-integral S harbor a gapless Bethe chain-like phase 
instead. The case of the TST is distinct from Haldane's 
analysis, however, owing to the different symmetry of the 
order. While the collinear Neel case is described by a vec- 
tor 0(3) NLaM, the spiral case is instead described by a 
NLaMwith a matrix SO (3) order parameter^. Here the 
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matrix may be constructed from the local spin order, 

Si ~ m(hi cosq • + n 2 sinq • r^), (112) 

where rii and ri 2 specify the plane of the spiral, with 
rii n 2 = 0, q the spiral wavevector, and m the amplitude 
of the quasi- static moment. One can construct from this 
the SO (3) matrix 



= (n 1 |n 2 |n 3 ) : 



(113) 



with 113 = fix x n 2 . 

If on short space and time scales, spiral order is 
present, we expect that an appropriate effective NLcrM 
action is given by 



">NL(jM 



(114) 

Y J dxdr {^Tr [d r O T d r O] + vTy [d x O T d x O] }. 

Note that, for a quasi- Id system with circumference L y , 
the effective coupling constant g ~ c/L y <C 1 for large L y , 
with some constant two-dimensional coupling constant c. 

Famously, in Haldane's analysis of spin chains with 
a vector 0(3) order parameter, the naive NLcrM action 
must be supplemented by a topological term™. Topology 
of the order parameter is also important here, but its na- 
ture is rather distinct from Haldane's case. For clarity, we 
compare and contrast the two situations here. The vector 
0(3) order parameter comprises a manifold isomorphic to 
the sphere S 2 . Its topology is summarized by the homo- 
topy groups U^S 2 ) = and U 2 (S 2 ) = Z. The former 
implies that there are no non-trivial loops on the sphere, 
and correspondingly no singular point defects in two di- 
mensions. The latter, second homotopy group implies 
that there are classes of non-trivial smooth configurations 
of the order parameter in two dimensions, parametrized 
by an integer. These configurations are skyrmions, lack- 
ing any singularity. Because of the lack of any singu- 
larity, the skyrmions appear in a continuum limit of the 
0(3) vector NLcrM , and modify the physics of the NLcrM 
through a topological #-term, which gives a phase factor 
to configurations with non-zero skyrmion number. Based 
on this NLcrM with #-term, Haldane postulated distinctly 
different behavior for integer and half-integer spin chains. 

In the matrix SO (3) case, the order parameter man- 
ifold is S 3 /Z 2 , and the corresponding homotopy groups 
are IIi(S 3 /Z 2 ) = Z 2 and U 2 (S 3 /Z 2 ) = 0. The trivial 
second fundamental group means that non-singular con- 
figurations of the order parameter have no topological 
distinctions. This implies that a continuum limit exists 
in which there are no topological defects and there is no 
topological term. Instead, the non-vanishing first homo- 
topy group implies that there are singular point defects 
in two dimensions, with an Ising character. Note that 
in our theory, these are point defects in space-time, or 
instantons. Such defects are well-known in classical two- 
dimensional non-collinear magnets, and are known as Z 2 
vortices^. They do not appear in the continuum NLcrM 



, but are allowed in a lattice theory. Instead, the proper 
way to treat them is to embed the continuum theory in 
a larger one in which the defects appear as operator in- 
sertions, with some fugacity and selection rules. This sit- 
uation is familiar from the Kosterlitz-Thouless analysis 
of the classical XY model, in which the naive continuum 
theory is just the Gaussian spin- wave line, and the de- 
fects are point vortices which are treated as a kind of 
Coulomb gaiPP. It occurs also in the quantum analysis 
of 2+1 dimensional collinear antiferromagnets, where the 
singular defects are hedgehogs or monopoles. The sepa- 
ration of these defects and the continuum theory is the 
basis of the theory of deconfined quantum criticality^. 

With this understanding, we may first consider the 
SO (3) matrix NLcrM with out a ny Z 2 vortices, which is 
simply described by Eq. (112). There is no topologi- 



cal term. This SO(3) NLcrM is, like all NLcrM's in two 
dimensions for non-abelian groups, asymptotically free. 
Lacking any quantum phase factors, we expect simply 
that it develops a gap at a length scale £ ~ e 9o/g ^ e~ Ly , 
and that order parameter (hence spin) correlations decay 
exponentially beyond this scale. The gap itself behaves as 
A ~ v/£. Note the difference from Haldane's case, where 
the term, which is non-trivial for half-integer spin, fun- 
damentally alters the behavior of the continuum NLcrM, 
leading to gapless behavior in the half-integer spin case. 
Here, there is no topological term, and the system is al- 
ways gapped with exponential spin correlations. 

Now we can consider the role of the Z 2 vortex instan- 
tons. Such a vortex is described in the field theory by 
an operator, ip, which inserts the vortex at a particular 
space-time point. It is crucial to consider the quantum 
numbers of a Z 2 vortex, i.e. how the operator ip trans- 
forms under physical symmetries. The relevant opera- 
tions are time-reversal, translation, and inversion. It can 
be argued (we discuss this in Appendix[D| that the vortex 
operator is invariant under time-reversal and translations 
along and transforms under the other two operations, 
translation along x, T x , and inversion, P, according to 

T x : x -> x + 1, -> (-1) L ^, (US) 
P : x -> -x, y -y, ip (-l) Ly ^. (116) 

From the above properties, we see that for odd L y , ijj has 
the transformation properties of a staggered dimerization 
operator. In general, two operators with the same sym- 
metry are expected to have non-zero overlap in the op- 
erator sense, and their correlations will be proportional. 
Thus, for odd L y , the Z 2 vortex operator %jj can be viewed 
as a staggered dimerization order parameter. 

Let us consider the correlations of ip. Its two-point 
correlation function is obtained by inserting two Z 2 vor- 
tices in the system at separated space-time points. When 
they are widely separated, the result should be just the 
pro duct of two independent Z 2 vortices. Naively, using 
Eq. ( 112 ), such a vortex has an action which diverges log- 



arithmically with the system size. However, its effective 
action is expected to be finite, due to the vanishing order 
and stiffness beyond the scale £. Roughly, the effective 
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FIG. 23: Entanglement entropy at (a) R = 0, (b) R = 0.2 
and (c) R — 0.7 at zero field. The oscillatory behavior with 
periodicity 2 shows the dimerized ground state. 



action for a single vortex is thus obtained by replacing 
the system size by £, so S v ~ ^ 9o/g 2 - Then we 

expect that 



lim (?/>OM0)) - e~ 2 ^ - e -W<? 2 - e~ cL y 



(117) 



with some constant c. The saturation to a finite value as 
x — » oo implies 7^ 0, and hence, for odd L y , the exis- 
tence of staggered dimer order. For even L y , there is no 
connection of Z2 vortices to dimerization, so although the 
former are present, the system forms simply a featureless 
gapped state. 

We can probe into this state by measuring the entan- 
glement entropy in DMRG for a range of anisotropics at 
zero field. We show this in Fig. |23j where an oscillatory 
behavior of period 2 gives clear evidence of the dimerized 
phase described above. 



B. Gapless states in low but non-zero field 

As argued in the previous subsection, the ground state 
in zero field is a non-magnetic dimerized state with a 
gap to all excitations. As a consequence of the gap, the 
ground state is unchanged by application of a sufficiently 
small field. The ground state changes when the field is 
large enough that a state with non-zero spin crosses the 
energy of the spin zero ground state. Generally, if the 
transition to a state of non-zero magnetization occurs 
continuously, we can think that the state with non-zero 
magnetization consists of a dilute set of elementary exci- 
tations above the zero field ground state. 

We must consider therefore the elementary excitations 
of the dimerized state, and in particular those which 
carry non-zero spin (as these couple to the field). The 
most important such excitations are the topological soli- 
ton excitations which are characteristic of the broken 



Ising symmetry of the dimerized state. Such solitons 
are domain walls, connecting the two distinct dimerized 
ground states. As is well-known from the study of the 
Majumdar-Gosh chain 45 , solitons of this type carry spin, 
and in particular for the TST, one can readily argue that 
the solitons carry half-integer spin, namely S z = 1/2, 
3/2, as shown in Fig. 24 Both values of the spin are 



possible, and generally differ in energy. The solitons are 
topological excitations insofar as they are non-local: they 
cannot be created by the action of any local operator on 
a dimerized ground state. In addition to the topological 
soliton excitations, non-topological excitations carrying 
spin S z = 1 also exist. They can be visualized either by 
replacing a singlet dimer by a triplet of aligned spins, or 
as a bound pair of S z = 1/2 solitons. 

Generally, if the magnetized state is realized as a di- 
lute system of non-topological S z = 1 triplons, then the 
dimerization is not disrupted and must persist for M > 0. 
Numerically, however, the dimerization appears to be dis- 
rupted at all non-zero M . We will assume henceforth 
that the magnetized state (at small M > 0) should be 
regarded as a collection of topological soliton excitations, 
and neglect the S z = 1 triplons. 

In general, the excitations can be characterized by spa- 
tial quantum numbers in addition to spin. For an excita- 
tion localized in x in the TST, we may consider the trans- 
formations under translations along y, T y: and under in- 
version, P. From Fig. 24 , it is clear that the S z = 3/2 



soliton is invariant under both. However, this is not the 
case for the S z = 1/2 soliton, which has additional struc- 
ture. In general, out of the three non-dimerized spins in 
the "core" of the domain wall, we can form three linearly 
independent states with S z — 1/2, 



\m) = 



V3 




(118) 



where £ = e 27r / 3 and m = 0, ±1. These are simply mo- 
mentum eigenstates along the 3 site chain. The state |0) 
is invariant under the T y and P operations, while the chi- 
rality eigenstates |±) form a two-dimensional irreducible 
representation. In general, the chirality states would dif- 
fer in energy from the scalar one. If we crudely model the 
soliton core as a three-site antiferromagnetic Heisenberg 
chain, then we see that the chirality states have lower en- 
ergy, so we expect that the elementary solitons take this 
form. Consequently, there are two chirality "flavors" to 
the S z = 1/2 solitons. 

To understand the impact of the solitons, we will need 
the relation between the microscopic lattice operators 
and those which describe the solitons. The simplest to 
consider is the dimerization operator, or the bond kinetic 
energy, B x ^ y = S x ^ y • S x +i, y . This is negative on singlet 
bonds and has zero average on bonds with uncorrelated 
spins. In a ground state, it oscillates with period 2 in the 
x direction. However, the singlets are shifted over by one 
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FIG. 24: (Color online) Toy picture of solitons. The blue 
shade, covering thr ee si tes, in (a) corresponds to the S z = 
1/2 soliton in Eq. (118) while the three blue arrows in (b) 
corresponds to a single S z = 3/2 soliton. 



sublattice on crossing a soliton, so 



B XtV ~B + (-i) 



x+N(x) 



(119) 



where B is the non-zero average, and eo is the ampli- 
tude of the bond modulation. We have defined N(x) = 

T, x > <x ai, x > a +,x> + a -,x> a -,x> + a t,x> a 3,x» Which 18 the 

number of solitons to the left of the position x. The 
N(x) factor accounts for the shift in the singlet position 
on crossing each domain wall. 

Next we turn to the spin density operator S xy . We 
are interested in its action on states which consist of a 
low density of solitons. It is helpful to consider a carica- 
ture of these states in which solitons are described by a 
wavefunction which is a product of columns of singlets, 
spaced by occasional non-singlet columns with either the 
chiral S z = 1/2 form, or fully aligned S z =3/2 spins, as 
shown in Fig. 
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If the operator S z y acts on a column 
x which is part of a singlet, it converts that singlet to an 
S z = triplet state. This triplet costs a non-zero energy 
equal to the zero field spin gap, and having S z = gains 
no energy back from the magnetic field. Thus if we re- 
strict our description to a low energy one, below the zero 
field spin gap, we can simply take S z y to annihilate the 
state in this case. If, however, x is located at the posi- 
tion of a soliton, then S xy gives back a low energy state, 
which consists either of the original soliton or one with 
reversed chirality. Notably, in moving down the Id sys- 
tem, solitons alternate between odd and even columns of 
the lattice. Thus a non-zero spin is only measured when 
y acts on an even or odd site, if the number of solitons 
to the left of the position x is fixed. This lets us write 
the following expression for the spin operator, 



x,y 



c+AT(a 



'+,x 



+C Va +,x a -,x + C y a-,x a +,x + a t,x a s ,x] > ( 12 °) 

where 3j , (X — x are annihilation operators for chiral 
S z = 1/2 solitons, and a^^ x is an annihilation operator 
for an S z = 3/2 soliton. 

Finally we consider the spin raising operator, S+ y , con- 
taining the XY components of the spin. Acting on a site 



which is part of a singlet bond, the raising operator con- 
verts the singlet to an S z = 1 triplet, with amplitude 
=Fl/\/2 depending upon whether the site is the left or 
right member of that singlet. The triplet with S z = 1 
has overlap with the state of two adjacent S z = 1/2 soli- 
tons (as well as other states not in the low energy sector) . 
Simple algebra shows that, for instance 




(121) 



->|- 



(l+>|-> + !-)!+» 



where on the left hand side, \s) represents the singlet 
state, and the columns represent the three columns in 
the TST. One the right hand side, the state has been de- 
composed into soliton states, and the ellipses represent 
higher energy states. Here we took the triplet to reside 
in the middle row. The other triplets can be obtained by 
translation, as the chirality states are translational eigen- 
states. From th is co nstruction, we obtain the analogous 
relation to Eq. (|120|), 



St v ~ (-!)•+"<"> £ [C my < x al !X+i _ ir+ 



JV(as) 



m=± 



(122) 



The low energy excited eigenstates will not consist of 
localized quasiparticles but delocalized ones, as solitons 
may hop between columns of the same sublattice, i.e. 
even or odd x. As a consequence, the states are eigen- 
states of the x- momentum k x , which is defined modulo tt 
rather than the usual 27r, due to the doubled background 
unit cell of the dimerization. In the dilute limit we should 
consider only the states near the minimum energy of the 
corresponding energy bands. For the S z = 3/2 solitons, 
which are inversion symmetric, if this minimum is non- 
degenerate it must occur at h x = or k x = tt/2. We 
expect it to occur at the latter, k x = tt/2 value, owing 
to the dominant antiferromagnetic spin correlations. For 
the S z = 1/2 solitons, inversion symmetry implies in- 
stead that if the positive chirality (q = +1) soliton has 
minimum energy at k x = qo, then the negative chirality 
soliton has its minimum energy at k x = —qo- We are 
not aware of a general argument to fix the momentum 
<2o, however, and expect it is generically non-zero. We 
have checked this by a crude and uncontrolled variational 
calculation of the soliton dispersion, which indeed gives 
minimum energy states with opposite non-zero momenta 
for opposite chirality (this calculation gives qo — 7r/6, 
but we do not expect this to be accurate). 

With this in mind, we focus only on the minimum en- 
ergy states and take a continuum limit, writing 



±iqox 



a 3 ,x - ^(x)e^ x , 



(123) 
(124) 
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where i/j m (x) and *fr(x) are take n as s lowly varying con- 
tinuum boson fields. Then Eqs. (|12Q ) ( 122 ) become 



Si,y ~ [l + (-1)-+^)] [ J2 tit* 

m=± 

+ Y,^ m{2q ° X+ ^ V) ^U-m + ^^l (125) 
m 

Sl y ~ 2isi n(/ oE eim(29oa:+4fy)m (^) 2 

m 

+2co S(2 o(-l) x+JV(x) ^e im(29oX+ *^(^) 2 

m 

+2cosg (-l)* +N( *ty+^I. (126) 



We are now in a position to write down an effective continuum theory to describe the low magnetization state in 
terms of bosonic field operators ip m for S z = 1/2 solitons with chirality m and \I/ m for the S z =3/2 solitons, all taken 
near their band minima. By symmetry, it takes the form 

H low = Jdx{ A("^^ + ei/2-V2)V' m + * t (-^^ + e3/2-3V2)* + V[^ +) ^_,* t *]}. 

(127) 



Here V is a general potential of quartic order and higher 
in the fields, representing interactions of the solitons. We 
have dropped terms above which mix the different soliton 
species, e.g. ones which might annihilate one S z = 3/2 
soliton while creating 3 S z = 1/2 solitons. Most such 
terms, at least at low order, are prohibited by various 
symmetries, such as translation and inversion symmetry, 
at least for a generic incommensurate wavevector go for 
the S z = 1/2 solitons. 

Consider increasing the magnetic field h from zero. 
The ground state remains the soliton vacuum, i.e. the 
dimerized state, until the energy of a state with non-zero 
solitons crosses the energy of the vacuum. Assuming re- 
pulsive interactions between solitons, this occurs when 
the energy of a single soliton vanishes, and this type of 
soliton will enter the system. We must compare the en- 
ergies ei/2 — h/2 and e 3 / 2 — 3/i/2, and see which vanishes 
first on increasing h. If the S z = 3/2 soliton energy is 
large, e 3 / 2 > 3ei/2, then the S z = 1/2 solitons will ap- 
pear, at h = 2ei/2- Conversely, if e 3 / 2 < 3ei/ 2 , then the 
S z = 3/2 solitons will appear, at h = 2e 3 / 2 /3. The crit- 
ical ratio e 3 / 2 Ai/2 = 3 is valid at infinitesimal soliton 
density, i.e. M — » + . At larger magnetization, interac- 
tions amongst solitons may become important, and will 
probably tend to disfavor the S z = 1/2 solitons further, 
since these must occur at a higher density and hence in- 
teract more strongly. Since in any case we do not know 
the energies e 3 / 2 , e^, we cannot actually use this criteria 
quantitatively. Instead, we simply consider both types of 
soliton liquids as possibilities, and determine their prop- 
erties at a phenomenological level. 



Let us consider first the S z = 1/2 case. Then we can 
neglect the ^ particle, which has an energy gap even 
when the ifj q solitons enter the system. The structure 
of the solitonic sta te is determined to a degree by the 
potential V in Eq. (127). By symmetry, it has the form 



F[n + ,n_,0] = -(n 



_) + 6n + n_ 



(128) 



With a > for stability, the state depends upon the coef- 
ficient b. If we assume b < a, then it is favorable for both 
solitons to enter the system in equal amounts, and the 
system forms a one dimensional Bose liquid of particles 
with two flavors. Owing to the strong quantum fluctu- 
ations in one dimension, this is a Luttinger liquid phase 
with two independent massless bosonic modes, associated 
to the two conserved densities. In the CFT terminology, 
this is a state with central charge c = 2. If instead b > a, 
it is preferable for the system to choose one state of soli- 
ton only. In this case there is a spontaneously broken 
discrete symmetry (inversion P), and only a single mass- 
less bosonic mode, or c = 1. We focus on the former case, 
b < a, which we argue describes the same phase as the 
semi-classical incommensurate planar state. 

To see this, we show that the spin correlations in the 
two- flavor S z = 1/2 soliton liquid have the same form 
as those in the Id incommensurate planar phase, de- 
scribed in Sec. |III E| In the soliton liquid, we can use 
the usual bosonization of bosons for each of the two 
species, ~ y/n^/2e~ l6rn , ~ n s /2 + d x <t> m /ir 

(and \E^\£ = 0), where (j) m is the dual field to the boson 
phase m . With this, we may conveniently represent the 
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non-local operator N(x) = n s x + J2 m ^m/ 71 "? where n s is 
the mean soliton density. Note since each soliton carries 
S z = 1/2 spread over the TST of width 3, the average 
magnetization per site is M = \n s /2 = n s /6. Then 

B x , y ~ S + e cos[(7r + 2(J)a; + ^], (129) 
S*„ ~ (l + cos[(7r + 25) a; + ^])(M+-^ 

5 [0 + -0_+2q o x+—y}J, (130) 
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+2 cos go cos[(7r + 25)x + V 9 ] ^e 1 



■)■ 



(131) 



Here 2 J = 7rn s = 27rM/3 and = 0+ + 0_. We can 
compare the above to the semi-classical result. In the 
semi-classical limit, the bosonic phases 0± are weakly 
fluctuating, while <p± and hence ip are strongly fluctu- 
ating. Then the dominant terms in the spin operators, 
with smallest scaling dimension, are those which do not 
contain any of the strongly fluctuating phases, 



2tt 

S z x , y ~ M + n s cos[0 + 2q x + —y], 



(132) 



2tt 



. -' -4 sin g e" sin[# + 2q x + —y], (133) 



where we defined = 6+ + 6- and 6 = 6+ — 0_. This can 



be directly compared to Eqs. (20) of Sec. HIE We see 



that the form of the spin operators is identical to that in 
the incommensurate coplanar state. Thus we can regard 
the S z = 1/2 chiral soliton liquid as another limit of the 
same phase. 

Let us turn to the case of the S z = 3/2 soliton liquid. 
As there is no chirality quantum number in this case, the 
state can be simply viewed as a Luttinger liquid without 
spin, and is expected to be described by a c = 1 theory 
of a single massless boson. We argue that this S z = 3/2 
soliton liquid is in fact another SDW phase very similar to 
the one obtained by the quasi-one-dimensional approach 
of Sec. |V C| While one might have expected to find the 
identical SDW phase in this way, we instead find that the 
S z = 3/2 soliton liquid is an SDW state with a different 
SDW wavevector, in particular with Q y = 0, contrasting 
with the value Q y = 27r/3 obtain from the quasi-ld ap- 
proach. If the S z — 3/2 liquid indeed occurs, therefore, 
we presumably require a phase transition to the other 
SDW state upon increasing magnetization. 

To observe the SDW structure of the S z = 3/2 soli- 
ton liquid, we again consider the spin correlations. Now 
we have no chiral solitons, i^rn^m = 0- This immediately 
implies that there are no low energy excitations with spin 
S z = 1 and hence no low energy content to the S ± oper- 
ators. Thus XY correlations decay exponentially in this 
phase, exactly as in the the SDW phase. To examine the 



S z correlations, we can bosonize the non-chiral bosons. 
This gives ^ ~ ^frT s e id , ~ n s + d x (p/n, with dual 

phases Now N(x) — n s x + y?/7r, and we note the 

relation between the magnetization and soliton density 
is changed to M = n s /2, since the solitons have spin 
S z = 3/2. We see then that 



(1 + cos[(tt + 2S)x + <p])(M - 



2tt J 



(134) 



Higher harmonics of the above cosine also appear in a 
more careful treatment. Note that the incommensurabil- 
ity is different in this case: 25 — wa s — 2ttM. Eq. (134| 
can be compared to the corresponding formula, Eq. ([82 



for the quasi-ld SDW state in the TST. We see that it is 
identical, save for the presence of a factor 2i\yjZ inside 
the cosine in the quasi-ld case. This shows that the two 
states have the same structure, save for a difference in 
the SDW wavevector, as mentioned above. 



VIII. DISCUSSION 

In this paper, we have presented a comprehensive anal- 
ysis of the field- anisot ropy phase diagram of the three- 
leg spin- 1/2 triangular spin tube, of interest primarily as 
an approximation to the corresponding two dimensional 
Heisenberg model on the anisotropic triangular lattice. 
Pronounced quantum effects, strongly deviating from the 
expectations based on classical analysis, occur through- 
out the phase diagram. In this section, we will discuss the 
implications of our results for two dimensions, and how 
robust these quantum effects are to other modifications 
to the model. 



A. Implications for two dimensions 

Throughout the paper we have commented on how re- 
sults obtained in the one-dimensional TST geometry ap- 
ply to the two-dimensional spin- 1/2 system. Here we 
summarize these connections, with particular attention 
to the phase diagram in 2d. With a few exceptions, 
the phases we obtained for the TST have straightfor- 
ward analogs in 2d, and consequently we expect the 2d 
diagram to be only slightly modified. 

For example, on the isotropic line, R = 0, away 
from very small field, all the phases we found are 
precisely those expected from the semiclassical analy- 
sis of Refs. 12131 We expect the semiclassical analysis 
to only work better in 2d, so the same coplanar and 
plateau states, and their incommensurate analogs for 
small anisotropy, < R <C 1, should occur there as 
well. We note that all coplanar states are intrinsically 
stabilized by quantum effects. 

Perhaps the most striking feature amongst these states 
is the 1/3 magnetization plateau, which extends well be- 
yond the semiclassical regime in our phase diagram for 
the TST. In fact, this plateau is further stabilized in 
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R=1-J'/J 



FIG. 25: Schematic phase diagrams for the two dimensional, 
S = 1/2 system. The shaded regions and boundaries con- 
taining full circles are based on preliminary DMRG results 
for circumference L y = 6,9 systems in addition to the TST 
with L y = 3. Other boundaries are drawn by hand using the 
considerations described in the text. Two possible schematics 
are drawn, differing in the extent of the cone phase. In (a) it 
is limited to the quasi- Id regime, while in (b) it extends to 
enclose the SDW state. The latter possibility is more clas- 
sical. Intermediate or more complex cases are also possible. 
See text for further explanations. The insets show plateau 
width Ah/ J as a function of R. 



2d, and we expect it to extend for all values of spatial 
anisotropy < R < 1. This expectation is based on 
both analytical arguments and additional DMRG stud- 
ies for cylinders of width L y = 6, 9 beyond the TST with 
L y = 3. Our results for the plateau width, Ah ~ 0.8J at 
R = 0, agree well with previous numerical studied 3 - ^ 33 * 46 *. 

Several experimental spin- 1/2 materials with the tri- 
angular lattice structure have indeed been observed 
to support a 1/3 magnetization plateau, including the 
well documented material Cs2CuBr j 47 * 48 * as well as 
Ba3CoSb20g, studied more recentl}@^. A notable ex- 
ception is CS2CUCI4, which is isostructural to Cs2CuBr4, 
but does not exhibit a magnetization platearP. In our 
opinion, as explained in detail in Ref. |4j the plateau is 
destabilized in this case by three-dimensional coupling, 
which is stronger (relative to the appropriate J) in the 
Cs-based magnet in comparison with the Br-based one^, 



with perhaps strong Dzyaloshinskii-Moriya (DM) inter- 
actions in CS2CUCI4 playing an additional role^l 

The SDW phase dominates a large fraction of the phase 
diagram for the TST. This is an entirely quantum phase 
(since it requires modulation of the length of the static 
moments), which in 2d exhibits incommensurate collinear 
long-range order along the field direction. Being of quan- 
tum origin, one may wonder whether the SDW persists 
into 2d. Based on renormalization group arguments, dis- 
cussed extensively in Ref. 4, we know that the SDW in- 
deed must exist in the quasi- Id regime, J' <C J, when 
inter-chain correlations are relatively weak. We expect 
that the region occupied by the SDW may be somewhat 
curtailed in 2d relative to that in the TST, but that it still 
is quite large. This is based on intuition and numerical 
evidence that inter-chain correlations remain suppressed 
for relatively large J' due to frustration. 

Experimental verification of this novel magnetic state 
is clearly called for. In this regard we would like to 
point out a recent series of experiments on quasi- Id spin- 
1/2 material LiCuV04. While much of the interest in 
this material stems from the high-field nematic phase 
predicted^ and observed^ to occur near the satura- 
tion field, several experimental studied^^ have found 
strong evidence in favor of an incommensurate longitu- 
dinal SDW phase in the intermediate range of magnetic 
fields. To understand this finding better, it is important 
to realize that the inter-chain exchange in this material is 
of zig-zag (triangular) type albeit of predominantly fer- 
romagnetic sigrPn The considerations of Section |v| make 
it clear that the SDW phase is not sensitive to the sign of 
inter-chain J' and should appear in the model with fer- 
romagnetic J' as well, see for example Ref. 60 for explicit 
calculations. We thus would like to posit that a recent 
neutron scattering study^, which observed longitudinal 
spin fluctuations but no transverse ones, is very much 
consistent with SDW phase scenario. Like the spin ne- 
matic phase, which is expected to occur at much higher 
magnetic field, the SDW phase too does not support low- 
energy transverse spin excitations. It would also be in- 
teresting to seek evidence of an SDW state in Cs2CuBr4. 

The above aspects of the TST and 2d phase diagrams 
are qualitatively similar. Qualitative differences are ex- 
pected at low and high fields. At zero field, the TS T ex- 
hibits a dimerized phase, which we attribute (Sec. VII) 
to quantum fluctuation effects specific to one dimension. 
In 2d, most of the zero field line should exhibit incom- 
mensurate spiral order, with a small region of collinear 
antiferromagnet at small J 7 / J, as argued in Ref. [5] At 
high field, near saturation, where the TST shows both 



coplanar and cone phases, we saw in Sec. IV D 2 that 
in 2d only the coplanar state occurs. This is a rather 
surprising result, since the coplanar state might be con- 
sidered more quantum than the cone. This observation 
poses a tricky problem of connecting the limit of field ap- 
proaching saturation at fixed small J', where the copla- 
nar state is expected, to the limit of vanishing J' at fixed 
field slightly below saturation, where we instead expect 
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a cone state. In 2d, therefore, a phase boundary must 
emanate from the saturation point at J' = 0, and we do 
not presently understand where this boundary extends 
to. 

Putting together all these considerations, we can con- 
struct schematic phase diagrams for two dimensions. The 
two simplest possibilities we could construct are shown 
in Fig. [25] The quasi-ld analysis, which was carried out 
directly in 2d in Ref. 4, demands the cone, SDW, and 
plateau phases at non-zero field and small J' /J. It also 
requires a collinear anti-ferromagnetic state at zero field 
and small J' / J . This collinear state is expected to be 
rapidly destroyed in favor of the SDW as the field is im- 
posed. It is likely to become canted as it does so, but in 
the absence of a detailed description of this narrow region 
descending from the collinear antiferromagnet at zero 
field, we label it "quasi-collinear" in the figures. Near the 
isotropic line, the semi-classical description requires com- 
mensurate (C) planar and incommensurate (IC) planar 
states, as well as the 1/3 plateau. Finally, near satura- 
tion, the dilute spin- flip approach becomes exact, and the 
solution of the BS equation required the IC planar phase. 
The shaded phases and the boundaries containing cir- 
cles are taken from preliminary DMRG results for more 
two-dimensional system with L y = 6, 9 lattice spacings 
around the circumference. The remaining phase bound- 
aries are drawn arbitrarily to connect the known regions 
demanded by the above reasoning in the simplest possi- 
ble manner consistent with scaling. The principle uncer- 
tainty in the diagrams is the extent of the cone phase. We 
expect it to occupy a relatively small portion of the phase 
space, despite the fact that it is the classical ground state 
everywhere below saturation except on the R = line! 



In the first schematic, Fig. 25 1, the cone state occupies 
the minimum possible area, while a more semi-classical 
situation might be as shown in Fig. [25b. 



1. Comparison to other work 

It is interesting to compare our results to those of Tay 
and Motrunich 33 , which is the only other comprehensive 
study of the full anisotropy-field phase diagram of which 
we are aware. We caution that a strict comparison is not 
possible because both their and our predictions for 2d are 
somewhat schematic, being based on conjectural extrap- 
olation of results for the Id TST (us) and finite clusters 
(them). Nevertheless, one notices immediately similari- 
ties between their schematic 2d phase diagram, Fig. 10 
of their paper, and our Fig. [3] First, the region near the 
isotropic line is in both cases quite close to semiclassical 
predictions. Small differences appear at low fields, where 
indeed quantum effects of the finite systems studied in 
both works are probably maximal. Second, near the sat- 
uration field, they also find a wide range of incommen- 
surate planar phase (called incommensurate V in their 
study). Our analytical BS analysis indicates that this 
phase in fact extends over the full range of anisotropy, 



a fact which was not resolved in their diagram. Third, 
both studies indicate the robustness of the 1/3 plateau. 
As already mentioned above, our results for the width 
of the plateau Ah ~ 0.8 J at the isotropic point R = 
agree well with those of Refs. l32|33l The more recent ex- 
act diagonalization stud}S^ predicts smaller width, about 
0.5J, but this is based on extrapolating Ah from small- 
size clusters. For R > 0, Ref. |33]is the only one we can 
compare with, and qualitative agreement is quite good. 
Our DMRG work completes the phase diagram, demon- 
strating the 1/3 plateau existence for all J' > 0. 

The major distinction between the two works is in our 
finding of the SDW state in a wide field anisotropy range, 
where Tay and Motrunich postulate separate spin liquid, 
spiral (corresponding to our cone state), and quasi-ld 
regimes. In our work, renormalization group arguments 
rather clearly establish the SDW phase in the small J' / J 
regime in 2d, which is the quasi-ld region of Tay and 
Motrunich. We think it likely that even in 2d, the SDW 
phase extends to R « 0.5. 



B. Suppressing the quantum effects 

As remarked above, we predict two types of quan- 
tum states - coplanar phases and collinear SDWs - in 
the 2d S=l/2 model. While remarkably robust in this 
case, these quantum phases can be suppressed by other 
changes to the model: larger spins 5 > 1/2, three- 
dimensional coupling, and Dzyaloshinskii-Moriya (DM) 
interactions. 



1. Higher spin 

We first consider 5 > 1/2, and find that the quan- 
tum phases are strongly suppressed. We begin with the 
vicinity of the saturation field. In Sec. |IVD~2} we showed 
that for 5 = 1/2 the system forms a coplanar state in 
this limit for all < J' / J < 1. This is surprising since 
except for the isotropic case, the coplanar phase is not 
a classical ground state. Using the calculations sketched 
below, we find that with increasing 5, the classical re- 
sults are recovered, with the coplanar phase restricted to 
increasingly narrow region near the isotropic limit, where 
it occurs due to classical degeneracy. 

To do so we use the representation below^ , which is 
more convenient than the Holstein-Primakoff one: 



s} = 

S' = 



2S[1 + 
S-b%. 



(K s - l)b%]b r , 



(135) 



where K s — yj\ — 1/(2*5). This expression reproduces 
the matrix elements of spin raising and lowering opera- 
tors between states with different magnetization exactly 
within the two-magnon (two spin flip) subspace. The ad- 
vantage of this form is t hat it requires no 1/5 expansion. 
Note that for 5 = 1/2 (p5]) reduces to p6l), thanks to 
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the hard-core condition (b r ) 2 = 0, while for large 5 ^> 1 
we recover Holstein-Primakov asymptote K s ~ —1/(45). 
Note that for 5 > 1 the hard-core constraint is not re- 
quired and as a result the [/-term is absent from the 
two-magnon HamiltoniarPH The Hamiltonian within the 
two-magnon subspace retains the form in (28) but now 



the interaction term is a bit more complicated, 

V(k, k', q) = l - ( J(q) + J(k + q - k')) (136) 
-SK S (j(k + q) + J(k' - q) + J(k) + J(k')) , 

J(k) = 2Jcos[kJ +4J / cos[y] cos[^^]. 



Numerical solution of the BS equation (44) for the two- 



dimensional triangular lattice, which proceeds along the 



same lines as in Sec. IV D finds that for higher spins 
5 > 1, near the saturation field the coplanar phase near 
the isotropic limit is limited to a region J' > J' CT > 0, with 
a cone phase obtaining instead for J' < J' cv . The critical 
value monotonically increases with 5, taking the values 
J' cv /J « 0.1,0.5,0.61 for 5 = 1, 3/2 and 2, respectively. 
These findings show that the absence of the cone state for 
5 = 1/2 found here is a very unusual feature of the most 
quantum case. Larger, more classical spins, do recover 
the classically expected state, although still in a limited 
range of J' /J. 

We next turn to the SDW phase. Since this state is 
rooted in the one-dimensional limit, we consider just the 
limit of weakly coupled chains, for 5 > 1/2, and in partic- 
ular 5 = 1. We find that the SDW is completely absent 
in this case. 

To see this, we consider a magnetic field above the 
lower critical field h& needed to overcome the non-zero 
Haldane gap (A a= i « 0.41 J for J' « J). This turns the 
gapped (and, essentially, decoupled - see Refs. l62|63l 
spin-1 chains into critical Luttinger liquids^ * 64 * 65 *. It 
turns out that these critical chains are characterized by 
a Luttinger parameter K = l/(47r7£ 2 ) > 1 for all val- 



ues of the magnetic field above the gap-closing /ia - 64 , 65 1 
This immediately implies that the scaling dimension of 
the longitudinal spin density operator S*_ 2 §(x) in (72) 
is K > 1 as well, which makes inter-chain SDW coupling 
in ( [T6] ) (which has twice this scaling dimension) strictly 
irrelevant. As a consequence the SDW phase does not oc- 
cur in the quasi- Id limit. Since this was its most stable 
regime in the 5 = 1/2 case, it may well be that the SDW 
phase is totally absent for 5 = 1! It would be interesting 
to check this in future simulations. 

What replaces the SDW? The large value of K implies 
an increased tendency to spin ordering transverse to the 
field direction, and indeed the twist term (4th term in 
(75)) is instead always relevant, leading to stabilization 
of the cone state. This result is supported by analytical 
and numerical studies of the spin 1 zigzag ladder. For 
example, Ref. 66 finds a finite vector chirality (that is, a 
cone state) for all values of the magnetization in the case 
of J i — J 2 spin-1 chain, along the J\ = J 2 line. 



Note that above, we found that the cone state was also 
stabilized for small J'/ J in the vicinity of saturation. It is 
likely then that the cone phase evolves smoothly between 
the Id limit J' / J = + and the approach to saturation at 
finite J' I J. Moreover, the presence of the cone state at 
small J' implies the asbence of any magnetization plateau 
in that regime. The predictions appear quite similat to 
those of the semiclassical analysis of Ref. [3j which sug- 
gests that the full phase diagram for 5 = 1 might be 
well described semiclassically. It is clear that in partic- 
ular the 1/3 plateau must terminate at some finite (and 
perhas not particularly small) value of the J' / J ratio in 
this case. 



2. Three dimensional coupling 

Another experimentally-relevant modification of the 
spin- 1/2 Hamiltonian is three dimensional coupling. We 
consider the simplest case of unfrustrated antiferromag- 
netic inter-plane exchange interaction J" between iden- 
tical triangular layers. Provided the three dimensional 
coupling is unfrustrated, we expect that the particular 
form is not too important. Such an interaction is ex- 
pected to make the spin system more classical and thus 
to promote the classical cone state over the coplanar one. 

Considering again the regime near saturation, one may 
readily solve the BS equation, appropriately modified to 
the three-dimensional situation. We indeed find that 
high-field co-planar configuration changes to the cone one 
for sufficiently large J" / J ratio. When the triangular lat- 
tice is isotropic, J' — J, this occurs for (J" / J) cr ~ 0.2, 
in agreement with the calculation in Ref. [I6j Not unex- 
pectedly, the critical J" becomes smaller for weaker inter- 
chain exchange J' . For example, for J' / J = 0.75, as per- 
haps appropriate for Cs2CuBr4, we find (J" / J) cv ~ 0.15 
while for J'/J = 0.34 (the Cs 2 CuCl 4 case), (J"/J) cr « 
0.034. On e-dim ensional scaling arguments, described in 
Appendix B3 1 suggest that (J" / J) cv ~ (J' / J) 2 when 
J 1 1 J <C 1, in agreement with the numerical values listed 
above. 

In the Id limit, J' / J <C 1, introduction of unfrustrated 
J" I J <C 1 disfavors SDW order in favor of a cone phase. 
This is discussed in detail in Sec.V of Ref|4j Thus three- 
dimensional coupling, if unfrustrated, tends to remove all 
quantum features of the phase diagram. 



3. Dzyaloshinskii-Moriya interactions 

A variety of DM interactions can be present in anis- 
otropic triangular lattice systems, depending upon the 
crystal symmetry and microscopic details. This can lead 
to diverse effects which are difficult to discuss without 
being more specific. For the materials CS2CUCI4 and 
Cs2CuBr4, the symmetry allowed DM interactions were 
obtained and discussed in detail in Refj4| Here we de- 
scribe only the effects of the dominant DM term in those 
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materials, which can be written as 

^dm = • S x ^ y x (S x _i 52y+ i — S^+i) , (137) 

x,y 

in the notation of this paper, with the DM-vector D = 
Da. oriented along the crystallographic a axis, normal to 
the triangular planes. 

Though small, a non-zero D has significant effects in 
both zero field and when a magnetic field is applied nor- 
mal to the triangular plane, i.e. parallel to the DM- 
vector. In these situations, unlike the J' interchain cou- 
pling, it is not frustrated either by the dominant chain 
interactions J or by the applied magnetic field. It tends 
to favor the cone state (or a spiral in zero field), and can 
obliterate the more quantum coplanar and SDW phases 
completely if sufficiently strong in this field orientation. 
Indeed, with this field orientation, an arbitrarily weak 
DM coupling inevitably forces the state in immediate 
proximity to the saturated state to be a cone phase, for 
all values of J' / J . This occurs because the DM coupling 
splits the degeneracy of the two minimum energy spin 
wave modes, already at the single spin wave level, making 
a two-component condensate impossible when the spin 
flip magnons are sufficiently dilute. 

We note, however, that when the magnetic field is ap- 
plied normal to the a axis, i.e. in the triangular plane, it 
itself frustrates the DM interaction. In this situation, the 
DM interaction is largely ineffective and has only minimal 
perturbative effects on the spin correlations. These field 
orientations are therefore optimal for observing quantum 
effects. 



C. Experimental implications and future directions 

Our study indicates that a number of "quantum" or- 
dered states may be found in S = 1/2 anisotropic trian- 
gular lattice systems. These states are not so exotic as 
quantum spin liquids, and are well characterized by their 
symmetries and associated order parameters. They are 
instead quantum in the weaker sense that they cannot 
be obtained in the classical limit. Most notably, we ob- 
tained a SDW state whose order involves (quasi-)periodic 
modulation of the length of the spin expectation value, 
along the field direction. We suggest this state occupies 
a wide swath of the field- anisot ropy phase diagram, pro- 
vided perturbations to our model are not too strong. 

The particular material Cs2CuBr4 appears a good can- 
didate for the observation of the SDW state, since three- 
dimensional coupling is known to be relatively weak 
there, and experiments have already identified the 1/3 
magnetization plateau. Direct observation of the SDW 
would consist of observing the incommensurate order- 
ing wavevector evolving monotonically with field, for 
fields above and below the plateau, and correlating this 
wavevector with the average magnetization. We ex- 
pect it to approximately follow the Id relation, q = 



7r(l — M/M 3 ), away from the plateau. Given its Id ori- 
gin, one might well also expect that the inelastic spectra 
retain Id features, such as spinon continua, in the SDW 
state and even in the plateau state above the gap. Of 
course, at low energy, in the vicinity of the SDW wavevec- 
tor, we expect the collective phason mode to dominate. 
There must therefore be significant rearrangement of the 
spectra on passing from low to high energy. A more de- 
tailed understanding of the spectral evolution with en- 
ergy, field, and anisotropy may make an interesting sub- 
ject for future study. 

In Cs2CuBr4, many additional features suggestive of 
phase transitions were identified above the 1/3 plateau 
in the magnetization process with an in-plane fieldP^Our 
study indicates that few such transitions should be ex- 
pected in the pure J — J' model. Likely addi tional DM 
interactions (beyond the one given in Eq. ( |137[ )) and per- 
haps further- neighbor couplings are at play. Study of 
their effects is a possible avenue for more research. 

More generally, the richness and surprisingly quantum 
nature of field- anisotropy phase diagram of the relatively 
weakly frustrated triangular lattice suggests that the be- 
havior on more frustrated lattices such as the kagome and 
pyrochlore may be even more interesting. The methods 
used here should be helpful in attacking these problems. 
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Appendix A: Sine-Gordon model and 
commensurate-incommensurate transitions 

In this appendix, we summarize the Commensurate- 
Incommensurate Transition (CIT) within the sine- 
Gordon model, which appears in multiple places through- 
out the manuscript. We consider the sine-Gordon action 
in d + 1 dimensions, with the form 

I /x 

-\cos[n($-qx)} j, (Al) 

where $ is the sine-Gordon field. We can write an alter- 
native expression in terms of the shifted field, $ = $ — qx, 



35 



so that 



M) 2 



+5d x d — A cos 



(A2) 



with 5 = p x q. In general, large S prefers an incommensu- 
rate state, where the field d is non-uniform and unpinned, 
while for small J, a commensurate phase occurs, where 
d is pinned to a fixed value by the cosine term. The 
detailed nature of the sine-Gordon model depends upon 
dimensionality, so we treat the d = 1 and d > 2 cases 
separately. 



1. 



> 2: mean-field transition 



For d > 2, the fluctuations of the phase field d are 
small even in the absence of the sine-Gordon term, i.e. 
for A = 0. This can be seen from the fact that, al- 
ready at the Gaussian level, the free boson propagator 
is non-divergent at small momentum for d > 2. This 
implies that the fluctuations of d are bounded, and one 
can therefore treat the entire problem by a saddle point 
approximation. Moreover, one can show that fluctuation 
effects are negligible in the (quantum) CIT for d > 2. 
More formally, Z) = <i+l = 2 + lis the upper critical 
dimension for the CIT. 

Therefore in this case we may proceed by simply mini- 
mizing the action in Eq. (A2). The minimum action con- 
figuration is independent of the d—1 coordinates normal 
to x and r. This gives 



(A3) 



5 where L± is the system width in the directions normal 
to x, and (3 is the length of the imaginary time integra- 
tion. The one-dimensional energy is then 



dx {^(dJf + SdJ- Acos(m9)}, (A4) 



where p = p x . Notice that 5 only appears as 
a boundary term, which means that the energy de- 
pends on 5 only through the winding number, N = 

U(x = L) - ${x = 0)) f-. Consider the case N = 0. 



Then, the solution is uniform, i.e. $ = 27rfc/n, with 
k = 0,1,2.... With N = 1, one obtains a well-known 
soliton solution of the sine- Gordon modeP^, which reads 



$(x) = — arctan \ e ±rM 
n 



(A5) 



where xq is the location of the center of the soliton. Note 
that the soliton has a width w ~ p/A, and energy E ~ 
\fp\. This gives a critical value, 



(A6) 



S c = 4y>V7r, 



such that, for 5 < £ c , domain wall solitons cost positive 
energy and so, are unfavorable, resulting in a commensu- 
rate wavevector. For S > 5 C1 it is favorable for solitons to 
be present, and the minimum energy configuration will 
be an array of solitons which characterizes an incommen- 
surate phase. 



Eq. ( A6 ) defines the location of the CIT phase bound- 
ary. We may also discuss its critical properties. On 
the commensurate side, no solitons are present, which 
implies the winding number N = precisely, and the 
ground state energy and field configuration are indepen- 
dent of 5. Thus, there is no visible critical behavior in 
the ground state (hence in equal time correlations) in 
the commensurate phase. On the incommensurate side, 
however, the minimum energy configuration of (f(x) de- 
pends upon S. It can be considered as an array of soli- 
tons, whose main characteristic is the spacing £ between 
solitons. This spacing is determined by the balance of 
the negative energy to introduce a soliton (which favors 
many solitons with a short spacing) and the repulsive en- 
ergy of interaction between solitons (which favors large 
spacing) . The repulsive interaction is exponentially small 
in the separation £ in units of the width w. Hence the 
energy of the array is 



Eid E ld 



■£/w 



(A7) 



where c is an unimportant constant, and L/£ is the total 
number of solitons. Minimizing this over £, one finds the 
critical behavior, to leading logarithmic accuracy, 



w In 



Sc 



S-5 C 



(A8) 



for < S — S c <C 5 C . The presence of the soliton array 
implies that the average gradient of the phase $ is non- 
zero, which defines the incommensurability wavevector 



q 



2tt 



1 



HM*-* c |/<yi 



e(5-5 c ). (A9) 



The incommensurability q in the incommensurate phase 
gives the shift of the ordering wavevector from its com- 
mensurate value. Other critical properties at the CIT in 
d > 2 are readily obtained from the results above. For 
example, the ground state energy density is simply the 
saddle point value of Eid, which scales as 



E 



S-Sc 



|m(£-<yi 



e(s-Sc). 



(A10) 



2. 



1: quantum fluctuations 



In the case d = 1, fluctuations of the phase field cannot 
be neglected. This can be anticipated from the Gaussian 
level result that, in the absence of a sine-Gordon term, 
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the free boson Green's function is logarithmically diver- 
gent at small momentum, signalling large fluctuations of 
Hence we must deal directly with the 1+1-dimensional 
action, 



Sg 



dxdrl -(d T d) z 



P 



(d x $) 2 + 5d x $- A cos mi? 



(All) 



Once again, 6 is the coefficient of a pure boundary term, 
which simply counts the number of solitons in the system. 
A finite density of solitons will be generated, provided 
the energy of a soliton for S = is compensated by this 
boundary energy, which equals 2it5/n. Thus we need the 
energy of a soliton at 5 = 0, i.e. in the pure quantum 
sine-Gordon model. 

We estimate this as follows. The scaling dimension of 
the cosine term, A n , is easily calculated, and is equal to 



(A12) 



The cosine is relevant when A n < 2, and irrelevant if 
A n > 2. When it is irrelevant, there is no pinning of the 
phase field at low energies. A state of this type is known 
as a "floating phase" , and because of the lack of pinning, 
the state becomes immediately incommensurate for any 
non-zero J, i.e. 5 C — 0, and there is no CIT. 

When the cosine is relevant, then when (5 = 0, the 
phase is pinned at low energies, and the energy of a soli- 
ton is non-zero. We need to estimate this energy to locate 
the value S c which defines the CIT. We do this by renor- 
malization group (RG) arguments. Renormalizing out to 
a length £, the cosine is reduced by fluctuations by an 
amount proportional to £ _A ™, so A e ^ ~ X^~ An . For a 
possible soliton of width £, the energy cost is of order 



(A13) 



The actual soliton size is determined by optimizing this 
over £, which gives 



and thus an energy cost for the soliton of order 



e s ~ A 2 



" A " in 2 ) 



(A14) 



(A15) 



This energy should equal the energy gain 2it5 c /n from 
the boundary term at the CIT, which gives 



(A16) 



Note that this approaches the mean-field result of the 
previous subsection when A n — >> 0, and becomes very 



suppressed when A n — >> 2~ (since we must assume A < p 
for consistency of the treatment). 

We now turn to the critical behavior, which in 1+1 di- 
mensions is a storied problem in critical phenomena. It is 
sometimes referred to as a Pokrovsky-Talapov transition, 
due to the solution by those authors.^ We recapitulate 
the essence of the argument. As in the mean- field case, 
for 5 < 5 C1 there are no solitons in the system, and the 
ground state energy is independent of 5, i.e. there is no 
sign of criticality in any static quantity. However, the ex- 
citation gap for creating a soliton vanishes linearly with 
S c — S. For < S — S c <C 5 C1 we expect a low density 
of solitons to be present in the system, again determined 
by the balance of the (negative) single soliton energy and 
the repulsive soliton-soliton interactions. 

We must, however, in this case treat the problem quan- 
tum mechanically. In particular, we must consider the 
effects of interactions properly in the low density limit. 
In this limit, the kinetic energy and momentum of in- 
dividual solitons is vanishingly small, and well-known 
results for low energy scattering apply. In particular, 
for short-range repulsively interacting particles in one di- 
mension, the probability of transmission vanishes in the 
low energy limit. Thus effectively, regardless of the mi- 
croscopic strength of the interaction, or of its short dis- 
tance structure, the solitons behave at low densities as 
though they were hard core particles, which cannot pass 
one another. To model this behavior, we can treat the 
solitons as fermions. Interactions at longer distances be- 
yond the local hard core are weak and unimportant, so 
the fermions are effectively free. 

The free fermion problem is trivially soluble, so we can 
easily obtain the critical behavior. When 5 > 5 C , we sim- 
ply fill the negative energy fermion states to form a Fermi 
sea. The sine-Gordon model has Lorentz invariance, so 
the dispersion of the solitons must be relativistic, hence 
the energy for a single soliton is 



sol 



y/e 2 8 +v 2 k 2 - —5, 



2tt , 
n 



(A17) 



where the velocity v = \J p/k, and e s = 2ii8 c /n. The 
Fermi momentum hp is determined by the condition 
^soi = 0. It will be small near the CIT, so we may 
expand the relativistic dispersion into its non-relativistic 
limit 



2-7T k 2 

E so] (k F ) = (S-6 c ) + -^ = 0, 

n zrn 



(A18) 



with m = e s /v 2 . This determines the Fermi momentum 



4irm 



(S - S c ) 



1/2 



(A19) 



The density of solitons is just hp/ft, as usual for spin-less 
fermions, so the incommensurability is thus 



2tt kr, 



2k, 



q= = — - ~ y/5 — 5 c . 

n 7r n 



(A20) 



37 



The square-root behavior is quite distinct from the log- 
arithmic one in d > 2. We may also easily obtain the 
behavior of the ground state energy density, as the total 
energy of the Fermi sea, 



E 
L 



2m n 



" F dk 
(5 - 5 c f/ 2 e(5 - 5 C ). 



5 C ) 



(A21) 



Many more results, e.g. for correlations in the incom- 
mensurate phase, can be readily obtained from the free 
fermion formulation, but we leave this to the reader to 
discover for themselves in the literature. 



Appendix B: Detailed calculations of BS 



In this appendix, we present our solutions to the Bethe- 
Salpeter (BS) equation in Eq. (44). This equation ap- 
plies only near saturation field, where the system can 
be modeled as dilute (hard core) bosons. We substitute 
our ansatz, Eq. ( [46] ) ? into the BS equation. With the 
constraint equation, Eq. ( [45] ) ? which enforces s = 1/2, 
we obtain a set of linear equations for the constants 
which can be written in a matrix form as 
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where we have defined 



Tij (/c, k ; fi) 



g e(k + q) + e(k' -q) + Sl 
T(q) = (l,cosq x ,smq x ,cosq y ,smq y , 
cos(q x - q y ), sm{q x - q y )), 



(B2) 



(B3) 



and Vt (x \h — h sat | . Although the s are integrals over 
simple trigonometric functions and other known quanti- 
ties, e.g. the dispersion, these integrals are divergent in 
both one and two dimensions and must be treated with 
care. It is possible, however, to analyze them asymptoti- 
cally. Once these integrals are evaluated, we can solve for 
the constants Ai to obtain T(q) from our ansatz, Eq. (46). 
In the next two subsections, we take the reader through 
our asymptotic analysis. 



1. Asymptotic behavior of for the 2d case 

In this section, we calculate the r^-'s for the 2d case. As 
aforementioned, we are interested in performing asymp- 
totic analysis in the limit Q — » 0, as the full inte- 
grals are too complicated to evaluate fully. We can 
partition the integrals into the the first two subleading 
terms, Bij ln(fi) + Qj, where the constants 5, C are in- 
dependent of Q. We can consider two cases: one with 
the same incoming momenta, i.e. the cone phase with 
ri = T(Q, Q, 0) = T(-Q, -Q, 0), and the other with dif- 
ferent incoming momenta, i.e. the coplanar phase with 
r 2 = T(Q,-Q,0) + T(Q,-Q,-2Q). Here, the wave 



vector Q minimizes the dispersion relation in Eq. (29), 



which can now be substituted into Eq. (B2). After some 
algebraic simplifications, we obtain 



Wo dqx L ' 



dq y 



a + b cos q y 



(B4) 



The exact forms of a, b will depend on whether the in- 
coming momenta are same or different. In this appendix, 
we will only present our results for i = j = 1, in which 



case, we can integrate analytically over q y in Eq. (B4), 
and obtain the following 



Til 



p2tt 

/ dq x 
Jo 



1 



b 2 



(B5) 



To proceed further, we need to specify the exact form of 
a and b. 

1. Same incoming momenta: For the same incoming 
momenta, a, b take on the following form 

a = SI + J(2 + j' 2 - (2 - j' 2 ) cos q x ) 

b = J(-2j' 2 cos(q x /2)), (B6) 

where we define j = J'/J. The integrand diverges 
near q x = like l/q x in the limit Q — » 0, and thus, 
integral is logarithmically divergent. After some 
analysis, the integral takes on the following form, 



1 



27rj'V4" 



: ln(fi) 



ln(2/(4-/ 2 )) 



(B7) 
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2. Different incoming momenta: For this case, a, b are 
as follows 

a = J(ft + 2+ j' 2 + (-2+ j /2 )cos(27rxi) 

+ j'y/^-j' 2 sin(27rxi) 
6 = -2Jj'(j' cos(ttxi) + V 4 - j' 2 sin(7rxi). (B8) 

The integrand now has two divergent points at q x = 
and q x = — arccos(l — j' 2 /2)/ir. Therefore, in 
comparison with the previous case, the logarithmic 
term doubles, and the integral takes on the form 

Tn i_ 1 n(n )+ 2h » , ' 4 -^» (B.) 



2. for the TST case 

In computing the r^-'s for the TST, we turn the two- 
dimensional integral in the previous section into a single 
integral over q x and a sum over q y . As one can imagine, 
the asymptotic behaviors differs in the TST from the 2d, 
in that, in the limit ft — >• 0, the integrals diverge as 1/ y/Q. 
Therefore, the two subleading terms of the integrals are 
Bij/y/U + Cij, where again, B,C are independent of Q. 
We present our results for i = j = 1 for the two cases of 
the same and differing incoming momenta. 

1. For the same incoming momenta, we obtain the 
following expression 



rii = — , : - + — , + 0(y/ti), (BIO) 

6 VP ~ f + WO 30)j' 2 + 24 v /(j / - + U' ~ 24/ + ^iff^ + 36 



where again, j = J'/J. ing momenta, in which case, the integral evaluates 

to 

2. We now compute t\\ for the case of differing incom- 

I 



■O(Vfi). (Bll) 



3. Weakly coupled chains limit 



In this appendix, we analytically check the results 
of Sec. |IV D 2| in the limit of weakly coupled chains, 
J' <C J. Recall that the calculation was done for a full 
two-dimensional lattice. Hereafter, we will use Cartesian 
coordinates, (x, y), for convenience. In this limit, we can 
express the spin flip operator as a continuous function 
of x, which is along the chain direction, while keeping 
the chain index y G 2 discrete. Then, from Eq. (71), 
we write this operator as ^ y (x) ~ S^ n (x), where its low 
energy theory is described by the following action 



Sid = X) / dxdr[^l(d r - ^d 2 - 

-t(*Jiftc* y+ i + h.c.) + u^l%% 
+ V *t*t +i$y+1$y J B ( B12) 



The spin-flip (magnon) mass, m = 1/J, follows from 
the quadratic dispersion of the magnon mode near mo- 
mentum 7r in a fully polarized chain. Additional inter- 
action terms describe the hard-core constraint (u term) 
as well as the transverse (t = J' xy /2) and longitudinal 
(v = 2J' Z ) parts of the interchain exchange interaction 
J' . Note that t-term contains a spatial derivative with 
respect to x, which reflects the frustration of the inter- 
chain exchange by the triangular geometry. In addition, 
this term contains a factor of i from the staggered factor 
( — \) x = e l7TX in Eq. (71), and from the fact that x takes 
half-integer values on odd chains (see Eq. Fig. [TJa), 
and Appendix D6 of Ref. HJ. 



We can analyze each term of Eq. (B12) through simple 
dimensional analysis, which will deem all these terms to 
be relevant under RG. Denoting the spatial scale along 
x as L, we can conclude that r ~ L 2 , ^ y scale as 1/y/L, 
while thee three interaction terms, t, u and v, scale as 
L. Hence, these are relevant interactions and must be 
included in our analysis of the low energy theory. 
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We can Fourier transform Eq. (B12) and write the 



Hamiltonian that corresponds to this action 

k 2 

k 



Hid = ^tf t k (^+2tk x cos[y- / i)'I' k 

k 

^ E ^(k,k',q)^ +q 4_ q ^ k ^ k(B 13) 



2N 

k,k',q 

Here V(k, k',q) = V(q) = 2u + 2vcos[q 2/ ]. Note that 
while the range of is not restricted, — oo < k x < oo, 
that of k y is limited by the lattice, — ir < k y < tt. This 
single particle dispersion contain two degenerate moment 
at Qi = (-2tm,0) and Q 2 = (2tra,7r). 

The single particle dispersion has two degenerate min- 
ima, at Qi = (— 2tra,0) and Q2 = (2tm,7r). We can now 
compute the renormalized couplings Ti,^ in a similar 
manner as the previous subsections. However, we alter 
our ansatz of the BS equation, Eq. ( [44] ) ? to take the form 
T(q) = Aq + AiCOs[q y ], because the odd contribution, 
oc sin[q 2y ], vanishes under the integral as the denomina- 
tor in Eq. (44) is even for all combinations of incoming 



and transferred momenta. 

Computing Ti, T 2 requires one to solve two linear equa- 
tions for Aq, Ai, which involve 2d integrals over functions 
with denominators like [k 2 /m + 16m£ 2 sin 2 [k 2y /2] +Q] (for 
ri) and [(k a ,-4mtsin 2 [k 2/ /2])/m + 4mt 2 sin 2 [k 2y ]+^] (for 
T2). We first evaluate these integrals analytically by sep- 
arating out the leading terms in \n[Q/(mt 2 )], then taking 
the limit u — >• 00 to, again, enforce the 8 = 1/2 constraint. 
The expressions are as follows 



87rf 

8irt 



i+jr 



(I + 57) lnT + 41n2 + 47(|ln2- 1) 
1 



lnT + 21n2' 



(B14) 
(B15) 



where 7 = v/(7rt) and T = 16mt 2 /Q. Given these forms, 
we can conclude that Ti > T 2 for 7 > j c = 3 In 2/(6 — 
4 In 2) « 0.644. Since we are considering the isotropic 
Heisenberg model, where 7 = 4/tt > 7 C , we observe that 
the coplanar fan state prevails over the cone state in the 
J'« J limit, in agreement with the full lattice approach 
in Eq. (48), once the parameters m,t,v are expressed in 



terms of exchange integrals. 

With this approach, we can also estimate the width of 
the planar fan state near saturation field through simple 
dimensional analysis of Eq. (B12). Since the chemical 
potential, \i = h sat — h, scales as L 2 and the t interac- 
tion scales as L, the phase boundary between the planar 
and the lower-field phase must scale as Ah ~ (J') 2 /J. 
This boundary separates the planar fan phase from the 
cone phase, a region in which a standard bosonization de- 
scription of Sec. [V] becomes appropriate. Details of this 
analysis are presented in Appendix |C 2| 

Similar reasoning allows one to estimate the stabil- 
ity of the planar fan state with respect to inter-layer 
coupling J", which is always present in real materi- 



als. It is clear that (non- frustrated) inter-layer cou- 
pling corresponds to adding a simple single particle hop- 
ping term between layers with a different z-coordinated 
J drdx^2 z J"(^ y z ^y 5 z+i + h.c.) term to the action in 
Eq. ( |B12 ). Such a term also scales as L 2 , which im- 
plies that the phase boundary between the planar and 
the cone phase in the J' — J" plane takes on a quadratic 
shape, J" - {J') 2 1 J. 



Appendix C: Additional one dimensional analysis 

The purpose of this appendix is to show that the TST 
geometry with 3-legs is unique in that the renormalized 
couplings generated through RG produce significantly 
different physics for N = 3 compared to that of N > 3, in 
the limit J' <C J. Moreover, we show that the arguments 
given below further support our claims in Sec. |VII| for the 
existence of a dimerized state near low field. Finally, we 
conclude this appendix with a more thorough analysis of 
the cone state near high fields. 



1. Zero field analysis by quasi- Id methods 

We start with the zero field case of Eq. ([T]) in the limit 
of decoupled chains J' <C J, where each Heisenberg chain 
can be bosonized using the Wess-Zumino-Witten SU(2)i 
theory, with central charge c = 1. In this theory, the spin 
operator can be decomposed into its uniform ~M. y (x) — 
Jr,v(x) +JL, y (x) and staggered N^(x) magnetizations 

S XtV a [M y (x) + (-l) x N y (x)} , (CI) 

and its scalar product can be written in the continuum 
limit 



-)■ (-l) x e y (x), 



(C2) 



where e y (x) is the staggered dimerization. With J' = 0, 
this theory describes the Luttinger liquid fixed point of 
the decoupled chains. The scaling dimensions of these 
continuum operators, M, N, and e, determine the rele- 
vance of each operator as it perturbs this fixed point. The 
uniform magnetization has scaling dimension 1, whereas 
both the staggered spin magnetization and the staggered 
dimerization have scaling dimension 1/2. These three 
continuum operators form a closed operator algebra with 
well-defined operato r pr oduct expansions (OPEs) used 
widely in liter at ure^^^f. For instance, the product of 
J# and N can be expanded as 



J a {x,r)N\x',r') 



ie abc N c (x,r) -i5 ab e(x,r) 



4tt (v(t — r r ) — i(x — x r ) + ao<r r ) ' 

(C3) 

where r is the imaginary time, v = irJao/2 is the spin 
velocity, and ao is the short-distance cutoff. 
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Let us now consider interchain Hamiltonian perturbing 
the decoupled Heisenberg chains, 

3 

V = J' E S v(*) (S„+i(ar) + S y+1 (a; - 1)) . (C4) 

y=l x 

Perturbation theory is formulated by expanding the par- 
tition function Z = J e ~ s o-S drV U p to quadratic order, 
i.e. 



Z~J e~ s ° 1 - jf V + jf jf V(ti)V(ts) 



(C5) 



with an implied short time cutoff a = ao/u. Here, T is 
the time-ordering operator. To utilize t his perturbation 
theory and the OPEs, we express Eq. (C4) in terms of 
continuum operators, Eqs. (CI) and ( |C2[ ), 



Vi = 



Vn 



y=l x 
3 

-atJ'J2J2 M yW -d x M y+1 (x), 



V* 



v 4 



(C6) 
(C7) 

3 i 

-^'EE N »W'o 9 X+i(4 (C9) 



y=l x 



y=l x 



where V = V\ + V2 + ^3 + V4. It is crucial to realize 
that the periodic boundary conditions enforced in the 7/- 
direction by the TST system, c.f. Fig. [2j allows us to 
rewrite any operator O as 



y=l x 



(CIO) 



= 1 X 



Using OPEs, one can show that the nearest neigh- 
bor chain couplings of the staggered magnetization and 
dimerization enter in the third power of J' , 



y=l x ^ ' 

(Cll) 

where J3 > and J3 ex (J f ) 3 . This is done by first 
generating d x N y -id x N y +i by quadratic in + terms. 
Next, this term is fused with V\ to generate the J3 oc 
(J') 3 interaction. The calculations are similar to those 
described in Refs. 4- 6|7Q|711 and refer the reader to these 
papers for more details. 

In a 2d system 4 5 , however, we find that the generated 
term is instead quartic in J', with interaction constant 
J4 ~ (J') 4 / J 3 and is of the opposite (negative) sign J 4 < 
in comparison with J3 above. It turns out that J3 ~ 
( J') 3 > is a feature of the N = 3 TST model only: 
wider tubes with N > 3 are anologous to the 2d case, 



where the renormalized couplings ~ (J f ) 4 /J 3 < 0. Note 
that this difference is important as it implies that spin 
tubes with N > 3 are not frustrated by the periodic BC 
along the ^/-direction. 

Going back to the N = 3 TST, both of the gener- 



ated interactions in Eq. (Cll ) are strongly relevant (scal- 
ing dimension 1) and scale to strong coupling under RG 
transformations. It would appear that beca use o f the 
greater numerical coefficient of € y € y +i in Eq. (Cll), it is 
the dimerized ground state that emerges from the com- 
petition in the strong coupling. However, this argument 
is not complete as it neglects the crucially important ef- 
fect of marginally irrelevant in-chain backscattering term, 
oc Jr • Jl, which in fact breaks the symmetry between 
the ~N y -N^+i and e y e y +i interactions in favor of the first 
one^. This outcome is not unexpected as it is well-known 
that in-chain marginal current-current interaction spoils 
the extended SU (2)r x SU (2) l symmetry of the Heisen- 
berg chain by subleading logarithmic corrections which 
modify chain spin correlations as follows^EHl 



(N y (x)N y (0)> = Qn[x]) 

(e y (x)e y (0)} = (ln[*])- 3 /V 



(C12) 



Essentially, the same mechanism promotes interchain 
N y • N y+ i interaction over that of staggered dimeriza- 
tions. In the infinite 2d lattice, this leads to the stabi- 
lization of the collinear antiferromagnetic phase 5 , which, 
however, is not possible in the TST geometry. 

It is important to realize at this point that the rele- 
vant J3 J2 y ~N y • N y+ i interaction, which describes non- 
frustrated coupling of staggered magnetizations on neigh- 
boring chains, changes the geometry of the system into 
that of a rectangular spin tube. The renormalized, rele- 
vant coupling, J3, become comparable to the intrachain 
exchange J under RG and forces Neel vectors Ni 5 2,3 to 
order into the familiar 120° pattern on every rung. Our 
Id reasoning stops at this scale, but further progress can 
be made by assuming that the spin tube with J3 ~ J 
can be accessed from the opposite limit of the strong 
rung exchange J± ^> fi^. In this limit, the spins on each 
rung form 3-spin triangles that interact via J± = J3, and 
are coupled to neighboring triangles by a weak exchange 
J. The ground state of each triangle is 4-fold degener- 
ate and is characterized by two quantum numbers, total 
spin 5 rung = 1/2 and chirality r, which is itself another 
pseudo-spin 1/2 object. The physical meaning of r is 
just a sense of either a clockwise or a counterclockwise 
rotation of the 'unpaired' spin- 1/2 in the ground state 
of the individual triangle. In other words, in addition 
to spin 1/2, the ground state now carries finite momen- 
tum ±27r/3 due to chirality. Focusing on this low-energy 
subset of triangle's states, one can derive spin-orbital 
Hamiltonian 75 



^s-o — jy ^ s rung 



(x) 



s rimg 



(x + l) X 



x[l + a N (r+T x+1 + r-r+ +1 )] (C13) 
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describing correlated dynamics of spins and chiralities. 
For the triangular ladder considered here, N = 3 and 
Of TV = 4. The presented arguments remain valid for any 
odd TV, however. See Ref. 75 for N = 5 and R ef. [76] for 
N > 5. A nalyti cal 1 ^ 1111 and numeric at- 75 1 78 1 79 1 studies of 
the model (C13) find dimerized ground state, in agree- 
ment with our consideration in Section |VII A| Fig. [23j 
which shows oscillatory behavior of the entanglement en- 
tropy for different values of R, represents clear evidence 
of the dimerized ground state. 

Finally, we conclude by discussing the way to generate 
an interaction of the uniform magnetizations from the 
next neighboring chains. This is done by fusing V\ in 
Eq. (C6) with itself, which yields, under Eq. (C5), 



SH 



MM 



(2J 1 



/\2 



2 



(M z y (x,r)M z y (x> >')> 

xm^m;, (cm) 



Because the result is converging, the integral of the y- 
th chain correlation function can be extended to the full 
x — vr plane. This, using important short-distance cut- 
off ~ aosign(r) and y = vr (see Ref. [71] for detailed 
discussion), leads to 



oo poo 

dx 

-oo 



dy(- 



1 



h.c.) 



47T. 



,(y + ix + a sign(y)) 2 

(C15) 

As a result we obtain for the amplitude of Sjmm = 
(J') 2 /(ttv), where v is magnetization dependent spin ve- 
locity. 



ence of an external magnetic field is given by 



V 

v cone 



^7cone ^2 



h.C. 



(C16) 



Even though the generated coupling constant is small, 
^7cone <^ J' / J <^ 1 5 this interaction does not involve 
spatial derivatives and has scaling dimension 2ir7Z 2 which 
approaches 1/2 as h — >• h sat . Thus, this is a strongly 
relevant term. 

In a 2d systenP^, 5j cone ~ (J') 4 / J 3 < as discussed 



in the previous subsection. (Note that (C16) is written 
in the 'sheared' system of coordinates.) When translated 
into Cartesian coordinates, it implies antiferromagnetic 
(positive) exchange interactions between spins on next- 
nearest chains at the same position x along the chairP. 
Crucially, as emphasized in the previous section, the TST 
geometry allows for a stronger renormalized coupling, of 
the order of Sj cone = h ~ (J') 3 /J 2 > 0. 

Th e difference is due to slightly different routes to 
(C16) in 2d and TV = 3 TST geometries. One can first 



show that, when you start from the original cone inter- 
action 



Kone — 7c« 



d^ ^t,y^xS^ y+1 



h.C, 



(C17) 



one can couple the derivatives d x S^ on the next-nearest 
chains y and y + 2, 



V" 



7< 



2 fY.j dx d ^y d * S ^y+2 + h.c. (C18) 



2. Cone state 

Now, turn on the magnetic field. When a large enough 
magnetic field is applied to the TST, the "twist" order, 
the fourth term in Eq. ( 75 ) , becomes more relevant than 



the SDW. This was discussed in previous papers for the 
two-chain ladder^ a s well as the 2d triangular lattice^. 
As both the SDW and the cone interaction amplitudes 
in (75) are of the order J', the relative importance of 



the two interactions can be estimated 5 from a compar- 
ison of their scaling dimensions, A saw = 1/(2tt71 2 ) and 
Acone = 1 + 2tt71 2 . These two dimensions are equal when 
2tt1Z 2 = (y/E—l)/2, which takes place at sufficiently high 
magnetization M « 0.6M sat . Because of rather steep de- 
pendence M(h) of the magnetization on the magnetic 
field near the saturation, this value of magnetization cor- 
responds to h « 0.9/i sat , see Fig. 2 in Ref. 5. A similar 
conclusion is obtained by comparing mean-field transi- 
tion temperatures of these two ordered states as functions 
of magnetization, see Ref. |4[ 

These arguments, however, are not complete be- 
cause they do not take into account the fluctuation- 
generated interactions between spin densities on next- 
nearest chains. The most important of these in the pres- 



This step parallels calculations leading to Eq. (Cll) with 



minor variation due to U(l) symmetry of the system in 
the presence of an external magnetic field. In this sit- 
uation the scaling dimension of the field is smaller 
than 1/2 which leads to a slightly different numerical pre- 
factor in the renormalization. However the functional 
dependence on J' remains the same. Secondly, for all 
N > 3 one also needs to generate 



^MM 



MM 



dx M;;M 



2/+2> 



(C19) 



which was des cribed in the end of the previous subsec- 

^ 2 /J. Fusing next 



tion, Sec. C 1 . Here, 5y 



7MM 



(C18) and (C19) together leads to the result (C16). In 



the N — 3 TST, however, the second step is not required 



due to (CIO), and we end up with a lar ger coupling of 
the order £ 7cone - (J') 3 / J 2 > in JC16| >- 

To compare the original V cone with the generated V^ one 
quantitatively, we can estimate the RG scale £ at which 
the coupling constant of the interaction beco mes o f the 
order one (in units of spin velocity v). For ( C17| ) this 
is, with logarithmic accuracy, £ C one ~ — ln (J / )/(2 — 
ln(J')/(l- 
2ttK 2 ) 



3) 



A cone ) 

-31n(J / )/(2 



2tt^ 2 ), while for (C16) it is £ 3 



We immediately conclude that 
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h < 4one for all values of 2nlZ 2 e (1,1 /2), i.e. that the 
generated cone interaction term ( C16| ) is more relevant 
than the bare one for all values of magnetization in the 
case of N — 3 TST. Similar consideration allows us to an- 
alyze the competition between the generated cone V^ one 
interaction and the SDW one, which is characterized by 
the RG scale 4dw ~ -2irlZ 2 ln(J' sin[<J])/(47rft 2 - 1). We 
find that 4dw < h for 1 > 2nlZ 2 > >/7 - 2 « 0.65, which 
corresponds to low-to-intermediate range of magnetiza- 
tion M > 0.25. At higher M, however, the modified cone 
interaction takes over the SDW one. (For the 2d case, 
the comparison is less conclusive as the result sensitively 



depends on numerical factors inside the argument of the 
logarithm 4 .) 

We now investigate the con seque nces of the strong 
J 3 = ^7cone interaction in Eq. (C16) for the TST prob- 



lem. In the high-field region where SDW fluctuations are 
suppressed, the Hamiltonian of the system is given by the 
sum of Hp in Eq. ( [69] ) ? the generated direct coupling V^ one 
in Eq. (C16), and the original cone interaction V cone in 
Eq. (C17), which now is a subleading one in comparison 
with (C16). With this, we perform abelian bosonization 



form of the interaction potential and arrive at the follow- 
ing expression, 



TST 



Js J dz{cos[/3(0i - 2 )] + cos[/3(0 2 - 3 )] + cos[/3(0 3 - 0i)]} (C20) 
J dxid^ + 2 ) sin[/3(0! - 2 )] + d x (0 2 + 3 ) sin[/3(# 2 - 3 )] + d x (0 3 + 0i) sin[/3(# 3 - 0i)]}. 



2 



r 



For J 3 ^> J 7 , which is the appropriate regime accord- 
ing to our RG arguments above, this potential is mini- 
mized by configurations with cos[/3(0 y — y +i)] = —1/2 
for all y. This allows for two different values of sine terms, 
sm[(3(0 y - Oy+i)] = ±y/S/2. In fact, different signs de- 
scribe states with different vector chiralities defined as 

K z y = (S y x S y+1 ) z ~ sm[/3(0 v - e v+1 )}. (C21) 

Thus, different signs of K z y correspond to different senses 
of rotation (clockwise or counterclockwise) of e z(36y as we 
go from one chain to the next. These chiralities also rep- 
resent useful order parameter describing two degenerate 
cone states^] 

To account f or th e subleading twist terms with spatial 
derivatives in (C20), we shift y — » y + vx, where v is 
determined by the requirement that in the new ground 
state, the bosonic field is twist-less, i.e. (d x 6 y ) = 0. 
Minimizing Hq + HjJ^ over v, we find 



v = -f3J f (sm[/3(6 y 



)]) 



(C22) 



This shows that the doubly-degenerate cone state is 
characterized by incommensurate transverse spin corre- 
lations, by virtue of the relation <S+ = (— l)^ e z/56>y — > 
exp[i(7r + v)x + i/36 y ]. Dep endin g on the spontaneously 
chosen vector chirality, Eq. (C21), transverse spin corre- 
lations are picked at either Qi iX = tt + v (for K z y > 0) or 
Q2,x — — tt + v (for Ky < 0) along the chain. 



Appendix D: Transformation properties of Z2 
vortices 



In this appendix, we address the transformation prop- 
erties of the Z 2 vortex instanton operator ip. We give sev- 



eral arguments. First, these properties have been implic- 
itly obtained in the case of a three leg spin tube, slightly 
different from the one studied here, in Ref.|8T| There, the 
authors explicitly evaluate the Berry phase contribution 
to the action for instantons on the lattice. Microscopi- 
cally, the instantons are associated with columns of spa- 
tial links along the x-direction of the cylinder (see below 
how this arises in another formulation). They showed 
that, due to the Berry phase, a single pair of instantons 
(an odd number of instantons cannot occur) is accompa- 
nied by a weight, 



e iS B p _ e 2KiS(x-x') 



(Dl) 



where x and x' are the locations of the instantons. For 
half-integer spins, this gives an oscillating factor equal to 
+1 or -1 if the separation between instantons is even or 
odd, respectively. From this we can extract the trans- 
formation properties. If we translate one of the instan- 
tons, x x + 1, we see that the weight in Eq. (Dl) 
changes sign. This requires ip —> —x/j, in agreement with 
Eq. ( 115 ). Under inversion, P, about a lattice site, the in- 



stantons, which live on the links, change from the even to 
odd sublattice of bonds and vice- versa. Inverting a single 
instanton, therefore, changes the pa rity of x, and hence 
also the sign of the weight in Eq. (Dl). Thus, again, 



ip is odd under inversion, in agreement with Eq. (116). 



Since the instantons do not move under time-reversal or 
translation along the invariance of ip under these op- 
erations is obvious. Thus, for the case L y = 3, for the 
model studied in Ref. |5TJ the symmetry of the instanton 
operator is determined as shown in the text. 

We turn now to an alternative derivation of the trans- 
formation laws, which gives the general result, and clar- 
ifies its generality. Here we follow the general strategy 
of Ref. |82j in which the Z 2 vortices are explicitly sepa- 
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rated from the smooth configurations of the SO (3) order 
parameter using a slave particle construction. This is 
achieved by writing the unit vectors defining the SO (3) 
matrix in terms of a "slave spinon" z a : 



ui + m 2 = e a f3Zf3(T ai z 1 , 



(D2) 



where the complex, two component vector z a is con- 
strained to have unit norm, ^ a z^z a = 1. This repre- 
sentation faithfully reproduces the orthonormality con- 
straints on the n^, but is two to one: the physical or- 
der parameter O is unchanged by the transformation 
z a — »• — z a . This is actually a gauge invar iance, since 
the transformation is made locally. The Z 2 vortex is a 
configuration in which, on encircling the center of the 
defect, z a returns not to itself but to — z a . 

As explained in Ref. |82j a low energy effective theory, 
appropriate to describe the regime with a local spiral or- 
der, as well as a quantum disordered phase, is a 2+1 
dimensional Z 2 gauge theory coupled to the spinon vari- 
ables z a . We refer the reader to Ref. [82] for details. The 
Z 2 vortex in this theory appears as a configuration of a 
spinon field which has a discontinuity z a — ^ — z a across a 
semi-infinite "cut" emanating from the vortex. This Z 2 
vortex is accompanied by an Ising vortex, the so-called 
"vison", which is itself a defect with a non-zero Ising 
gauge field crossing the same semi-infinite "cut" . In this 
way, the topological defects of the spiral magnet become 
identified with the visons of the Z 2 gauge theory. 

The discussion in the previous paragraph applies to Z 2 
vortices in two-dimensional space, which are particles in 
the 2+1-dimensional theory. We need to go from this to 
the description of instantons in the 1+1-dimensional the- 
ory obtained by applying periodic boundary conditions 
in the y direction. A 1+1-dimensional instanton can be 
viewed as an event in which a pair of Z 2 vortices is nucle- 
ated: one of them winds around the cylinder and finally 
arrives back at the other Z 2 vortex and annihilates it. 
We can, by the previous argument, consider the parti- 
cles nucleated and annihilated to be visons in the gauge 
theory. 



Such a process was considered in Ref. [83] (in the Sup- 
plementary Material), where it was shown that the oper- 
ator representing this process in the Ising gauge theory 
has the transformation properties in Eqs. (115), (116), 
i.e. this operator can be viewed as a staggered dimer- 
ization operator for odd L y . There, a rectangular lattice 
gauge theory was studied, but the basic physics is quite 
general. Let us consider the translation. We ask about 
the amplitude to first wind one vison around the cylinder 
at position x, and then wind another at position x + 1. 
The overall phase of the amplitude for both processes 
taken together gives the transformation property of the 
instanton operator under translation. The visons reside 
at the plaquette centers of the original lattice, and the 
winding trajectories form closed circles at fixed x, cir- 
cumnavigating the cylinder. Together, these two events 
form two such circles that enclose one column of sites 
in the lattice. The fundamental property of a vison is 
that it has a mutual statistical interaction with "electric" 
gauge charges, with the wavefunction acquiring a phase 
of 7r whenever one encircles the other. For a S = 1/2 
system, a unit gauge charge is present at every lattice 
site - this represents the physical spin at each site. The 
net effect of the two events together is that one vison is 
wound around each site of the lattice between the two 
circles, leading to an overall amplitude of (— l) Ly for the 
two processes together. Here, L y is the number of sites 
contain ed b etween the two circles. This gives the result 
in Eq. (115). Note that we may also roughly understand 



this phase factor by considering the smooth rotations of 
microscopic spins between the two contours, all of which 
rotate by 27r, and, due to their 5 = 1/2 spinor transfor- 
mation properties, each acquire a minus sign. A simi- 
lar argument shows that spatial inversion gives the same 
phase factor. Explicit calculations for these factors in the 
Ising gauge theory can be found in Ref. [83] Note that 
these arguments do not depend at all on the interactions 
in the model, just the presence of these symmetries and 
fundamental statistics of the particles. 
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